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Abstract. We outline the confluence of three novel theoretical fields in our modeling of Gamma-Ray Bursts (GRBs): 1) the 
ultrarelativistic regime of a shock front expanding with a Lorentz gamma factor ~ 300; 2) the quantum vacuum polarization 
process leading to an electron-positron plasma originating the shock front; and 3) the general relativistic process of energy 
extraction from a hlack hole originating the vacuum polarization process. There are two different classes of GRBs: the 
long GRBs and the short GRBs. We here address the issue of the long GRBs. The theoretical understanding of the long 
GRBs has led to the detailed description of their luminosities in fixed energy bands, of their spectral features and made 
also possible to probe the astrophysical scenario in which they originate. We are specially interested, in this report, to a 
subclass of long GRBs which appear to be accompanied by a supernova explosion. We are considering two specific examples: 
GRB980425/SN1998bw and GRB030329/SN2003dh. While these supemovae appear to have a standard energetics of 10^® 
ergs, the GRBs are highly variable and can have energetics 10^ - 10^ times larger than the ones of the supemovae. Moreover, 
many long GRBs occurs without the presence of a supernova. It is concluded that in no way a GRB can originate from a 
supernova. The precise theoretical understanding of the GRB luminosity we present evidence, in both these systems, the 
existence of an independent component in the X-ray emission, usually interpreted in the current literature as part of the GRB 
afterglow. This component has been observed by Chandra and XMM to have a strong decay on scale of months. We have 
named here these two sources respectively URCA-1 and URCA-2, in honor of the work that George Gamow and Mario 
Shoenberg did in 1939 in this town of Urea identifying the basic mechanism, the Urea processes, leading to the process of 
gravitational collapse and the formation of a neutron star and a supernova. The further hypothesis is considered to relate 
this X-ray source to a neutron star, newly born in the Supernova. This hypothesis should be submitted to further theoretical 
and observational investigation. Some theoretical developments to clarify the astrophysical origin of this new scenario are 
outlined. We turn then to the theoretical developments in the short GRBs: we first report some progress in the understanding 
the dynamical phase of collapse, the mass-energy formula and the extraction of blackholic energy which have been motivated 
by the analysis of the short GRBs. In this context progress has also been accomplished on establishing an absolute lower limit 
to the irreducible mass of the black hole as well as on some critical considerations about the relations of general relativity and 
the second law of thermodynamics. We recall how this last issue has been one of the most debated in theoretical physics in the 
past thirty years due to the work of Bekenstein and Hawking. Following these conceptual progresses we analyze the vacuum 
polarization process around an overcritical collapsing shell. We evidence the existence of a separatrix and a dyadosphere 
trapping surface in the dynamics of the electron-positron plasma generated during the process of gravitational collapse. We 
then analyze, using recent progress in the solution of the Vlasov-Boltzmann-Maxwell system, the oscillation regime in the 
created electron-positron plasma and their rapid convergence to a thermalized spectrum. We conclude by making precise 
predictions for the spectra, the energy fluxes and characteristic time-scales of the radiation for short-bursts. If the precise 
luminosity variation and spectral hardening of the radiation we have predicted will be confirmed by observations of short- 
bursts, these systems will play a major role as standard candles in cosmology. These considerations will also be relevant for 
the analysis of the long-bursts when the baryonic matter contribution will be taken into account. 


* We recall that “New perspectives in physics and astrophysics from the theoretical understanding of Gamma-Ray Bursts, I” was already published 

in COSMOLOGY AND GRAVITATION: X'* Brazilian School of Cosmology and Gravitation; 25''' Anniversary (1977-2002), M. Novello, S.E. 



INTRODUCTION 


In the last century the fundamental discoveries of nuclear physics have led to the understanding of the thermonuclear 
energy source of main sequence stars and explained the basic physical processes underlying the solar luminosity (see 
e.g. M. Schwarzschild [1]). 

The discovery of pulsars in 1968 (see Hewish et al. [2]) led to the first evidence for the existence of neutron stars, 
first described in terms of theoretical physics by George Gamow as far back as 1936 [3]. It became clear that the pulsed 
luminosity of pulsars, at times 10^ - 10^ larger than solar luminosity, was not related to nuclear burning and could be 
simply explained in term of the loss of rotational energy of a neutron star (Gold [4, 5]). For the first time it became so 
clear the possible relevance of strong gravitational fields in the energetics of an astrophysical system. 

The birth of X-ray astronomy thanks to Riccardo Giacconi and his group (see e.g. Giacconi and Ruffini [6]) led to 
a still different energy source, originating from the accretion of matter onto a star which has undergone a complete 
gravitational collapse process: a black hole (see e.g. Ruffini & Wheeler [7]). In this case, the energetics is dominated 
by the radiation emitted in the accretion process of matter around an already formed black hole. Luminosities up to 
10'*^ times the solar luminosity, much larger then the ones of pulsars, could be explained by the release of energy in 
matter accreting in the deep potential well of a black hole (Leach and Ruffini [8]). This allowed to probe for the first 
time the structure of circular orbits around a black hole computed by Ruffini and Wheeler (see e.g. Landau and Lifshitz 
[9]). This result was well illustrated by the theoretical interpretation of the observations of Cygnus-Xl, obtained by 
the Uhuru satellite and by the optical and radio telescopes on the ground (see Fig. 1). 

The discovery of gamma-ray bursts (GRBs) sign a further decisive progress. The GRBs give the first opportunity 
to probe and observe a yet different form of energy: the extractable energy of the black hole introduced in 1971 
(Christodoulou and Ruffini [11]), which we shall refer in the following as the blackholic energy^. The blackholic 
energy, expected to be emitted during the dynamical process of gravitational collapse leading to the formation of the 
black hole, generates X-ray luminosities 10^^ times larger than the solar luminosity, although lasting for a very short 
time. 

The extreme regimes of GRBs evidence new and unexplored regimes of theoretical physics. It is the aim of this talk 
to outline the progress achieved in understanding these astrophysical systems and the theoretically predicted regimes 
for the first time submitted to direct observational verification. 

It is a pleasure to present these results in Brazil. While sitting at the Casino de Urea, George Gamow and Mario 
Schoenberg in 1939 identified the basic process leading to the formation and cooling of a newly born neutron star (see 
). They called this process essentially related to the emission of neutrinos and antineutrinos the Urea process. It is a 
welcomed coincidence that, during the preparation of my talk at the tenth Marcel Grossmann Meeting [12], examining 
the data of the recently observed GRB 030329, we have received a confirmation of a scenario we have recently outlined 
in three papers giving the theoretical paradigms for the understanding of GRBs (Ruffini et al. [13, 14, 15]). 

We have clear evidence, first advanced in the system GRB980425/SN1998bw (Ruffini et al. [16], Fraschetti et al. 
[17]) and now confirmed also in the system GRB030329/SN1003dh, that there are in these systems three different 
components: 1) the GRB source, generated by the collapse to a black hole, 2) the supernova, generated by the collapse 
of an evolved star, 3) an additional X-ray source which is not related, unlike what is at times stated in the literature, 
to the GRB afterglow. In honor of the work done in the town of Urea by George Gamow and Mario Schoenberg, 
identifying in the neutrino emission of the Urea process the basic mechanism leading to the process of gravitational 
collapse and the formation of a relativistic compact star, we named these two X-ray sources URCA-1, the one formed 
in the system GRB980425/SN1998bw, and URCA-2, the one formed in the system GRB030329/SN2003dh. We shall 
now recall some of the main steps in reaching this understanding out of the GRB phenomenon and explore possible 
explanation of the origin of these two sources. 

We then turn to the analysis of the short GRBs. We first review some progress in the study of the general relativistic 
collapse of a shell of matter endowed with electromagnetic fields, which has been motivated by the study of the short 
GRBs. We then deduce from these theoretical developments some consequences for the interpretation of the mass- 
energy formula of the black hole, as well as some conceptual consequences for the relation between general relativity 
and thermodynamics. We also apply some recent progress on the solution of the Einstein-Maxwell-Blasov equations 
to the thermalization process occurring in the electron-positron plasma generated by the vacuum polarization process. 


Perez-Bergliaffa (eds.), AIP Conf. Proc., 668, 16 (2003), see Ref.[28]. 

^ This name is the English translation of the Italian words “energia buconerale”, introduced by lacopo Ruffini, December 2004, here quoted by his 
kind permission. 



The identification of the 
first black hole in our 
galaxy: Cygnus X-1 

• O = 1 (K' erg/s = 1 (H 7.0 

■ Absence uf pulsation due to 
uniqueness of Kerr-Newmann 
black holes 

• M > 3.2 V/ 0 

Leach & Ruffini, 1973 



Figure 1. Cygnus X-1 offered the possibility of identifying the first black hole in our galaxy (Leach and Ruffini [8]). The 
luminosity <I> of 10“^ solar luminosities points to the accretion process into a neutron star or a black hole as the energy source. 
The absence of pulsation is naturally explained either by a non-magnetized neutron star or a Kerr-Newmann black hole, which has 
necessarily to be axially symmetric. What identifies the black hole unambiguously is that the mass of Cygnus X-1, larger than 9Mq, 
exceeds the absolute upper limit of the neutron star mass, estimated at 3.2Mq by Rhoades and Ruffini [10]. 


We finally give a very specific theoretical prediction on the burst structure to be expected in short GRBs. We then 
proceed to the conclusions and some general considerations about this novel astrophysical scenario. 


THE ENERGETICS OF GAMMA-RAY BURSTS 

It is well known how GRBs were detected and studied for the first time using the Vela satellites, developed for military 
research to monitor the non-violation of the Limited Test Ban Treaty signed in 1963 (see e.g. Strong [18]). It was clear 
from the early data of these satellites, which were put at 150,000 miles from the surface of Earth, that the GRBs did 
not originate either on the Earth nor in the Solar System. 

The mystery of these sources became more profound as the observations of the BATSE instrument on board of the 
Compton Gamma-Ray Observatory (CGRO) satellite^ over 9 years proved the isotropy of these sources in the sky (See 
Eig. 2). In addition to these data, the CGRO satellite gave an unprecedented number of details on the GRB structure, 
on their spectral properties and time variabilities which became encoded in the fourth BATSE catalog [19] (see e.g. 


see http://cossc.gsfc.nasa.gov/batse/ 
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Figure 2. Position in the sky, in galactic coordinates, of 2000 GRB events seen by the CGRO satellite. Their isotropy is evident. 
Reproduced from BATSE web site by their courtesy. 


Fig. 3). Out of the analysis of these BATSE sources it soon became clear (see e.g. Kouveliotou et al. [20], Tavani [21]) 
the existence of two distinct families of sources: the long bursts, lasting more then one second and softer in spectra, 
and the short bursts (see Fig. 5), harder in spectra (see Fig. 4). We shall return shortly on this topic. 

The situation drastically changed with the discovery of the afterglow by the Italian-Dutch satellite BeppoSAX 
(Costa et al. [22]) and the possibility which led to the optical identification of the GRBs by the largest telescopes in 
the world, including the Hubble Space Telescope, the Keck Telescope in Hawaii and the VFT in Chile, and allowed as 
well the identification in the radio band of these sources. The outcome of this collaboration between complementary 
observational technique has led to the possibility of identifying in 1997 the distance of these sources from the Earth 
and their tremendous energy of the order up to 10^"*^ erg/second during the burst. It is interesting, as we will show in 
the following, that an energetics of this magnitude for the GRBs had previously been predicted out of first principles 
already in 1974 by Damour and Ruffini [23]. 

The resonance between the X- and gamma ray astronomy from the satellites and the optical and radio astronomy 
from the ground, had already marked in the seventies the great success and development of the astrophysics of binary 
X-ray sources (see e.g. Giacconi & Ruffini [6]). This resonance is re-proposed here for GRBs on a much larger scale. 
The use of much larger satellites, like Chandra and XMM-Newton, and dedicated space missions, like HETE-2 and, in 
the near future. Swift, and the very fortunate circumstance of the coming of age of the development of unprecedented 
optical technologies for the telescopes offers opportunities without precedence in the history of mankind. In parallel, 
the enormous scientific interest on the nature of GRB sources and the exploration, not only of new regimes, but also 
of totally novel conceptual physical process of the blackholic energy, make the knowledge of GRBs an authentic new 
frontier in the scientific knowledge. 









Figure 3. Some GRB light curves observed by the BATSE instrument on board of the CGRO satellite. 


THE COMPLEXITY AN SELF-CONSISTENCY OF GRB MODELING 

The study of GRBs is very likely “the” most extensive computational and theoretical investigation ever done in physics 
and astrophysics. There are at least three different fields of research which underlie the foundation of the theoretical 
understanding of GRBs. All three, for different reasons, are very difficult. 

The first field of research is the field of special relativity. As I always mention to my students in the course 
of theoretical physics, this field is paradoxically very difficult since it is extremely simple. In approaching special 
relativistic phenomena the extremely simple and clear procedures expressed by Einstein in his 1905 classic paper 
[24] are often ignored. Einstein makes use in his work of very few physical assumptions, an almost elementary 
mathematical framework and gives constant attention to a proper operational definition of all observable quantities. 
Those who work on GRBs use at times very intricate, complex and often wrong theoretical approaches lacking the 
necessary self-consistency. This is well demonstrated in the current literature on GRBs. 
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Figure 4. The energy fluence-averaged hardness ratio for short (T < 1 s) and long (T > 1 s) GRBs are represented. Reproduced, 
by his kind permission, from Tavani [21] where the details are given. 


The second field of research essential for understanding the energetics of GRBs deals with quantum electrodynamics 
and the relativistic process of pair creation in overcritical electromagnetic fields. This topic is also very difficult but for 
a quite different conceptual reason: the process of pair creation, expressed in the classic works of Heisenberg-Euler- 
Schwinger [25, 26] later developed by many others, is based on a very powerful theoretical framework but has never 
been verified by experimental data. The quest for creating electron-positron pairs by vacuum polarization processes in 
heavy ion collisions or in lasers has not yet been successfully achieved in Earth-bound experiments (see e.g. Ruffini, 
Vitagliano, Xue [27]). As we will show here, there is the tantalizing possibility of observing this phenomenon, for the 
first time, in the astrophysical setting of GRBs on a more grandiose scale. 

There is a third field which is essential for the understanding of the GRB phenomenon: general relativity. In this 
case, contrary to the case of special relativity, the field is indeed very difficult, since it is very difficult both from 

















































































































Figure 5. On the upper right part of the figure are plotted the number of the observed GRBs as a function of their duration. The 
bimodal distribution corresponding respectively to the short bursts, upper left figure, and the long bursts, middle figure, is quite 
evident. The structure of the long bursts as part of the afterglow phenomena of GRBs is illustrated in section . 


a conceptual, technical and mathematical point of view. The physical assumptions are indeed complex. The entire 
concept of geometrization of physics needs a new conceptual approach to the field. The mathematical complexity of 
the pseudo-Riemannian geometry contrasts now with the simple structure of the pseudo-Euclidean Minkowski space. 
The operational definition of the observable quantities has to take into account the intrinsic geometrical properties and 
also the cosmological settings of the source. With GRBs we have the possibility to follow, from a safe position in an 
asymptotically flat space at large distance, the formation of the horizon of a black hole with all the associated relativistic 
phenomena of light bending and time dilatation. Most important, as we will show in details in this presentation, general 
relativity in connection with quantum phenomena offers, with the blackholic energy, the explanation of the tremendous 
GRB energy sources. 

For these reasons GRBs offer an authentic new frontier in the field of physics and astrophysics. It is appropriate 
to mention some of the goals of such a new frontier in the above three fields. We recall in the special relativity field, 
for the first time, we observe phenomena occurring at Lorentz gamma factors of approximately 300. In the field of 
relativistic quantum electro-dynamics we see for the first time the interchange between classical fields and the created 
quantum matter-antimatter pairs. In the field of general relativity also for the first time we can test the blackholic 
energy which is the basic energetic physical variable underlying the entire GRB phenomenon. 

The most appealing aspect of this work is that, if indeed these three different fields are treated and approached 
with the necessary technical and scientific maturity, the model which results has a very large redundancy built-in. The 
approach requires an unprecedented level of self-consistency. Any departures from the correct theoretical treatment 
in this very complex system lead to exponential departures from the correct solution and from the correct fit of the 






















observations. 

It is so that, as the model is being properly developed and verified, its solution will have existence and uniqueness. 


GRBs and special relativity 


The ongoing dialogue between our work and the one of the workers on GRBs, rests still on some elementary 
considerations presented by Einstein in his classic article of 1905 [24]. These considerations are quite general and 
even precede Einstein’s derivation, out of first principles, of the Lorentz transformations. We recall here Einstein’s 
words: “We might, of course, content ourselves with time values determined by an observer stationed together with 
the watch at the origin of the co-ordinates, and co-ordinating the corresponding positions of the hands with light 
signals, given out by every event to be timed, and reaching him through empty space. But this co-ordination has the 
disadvantage that it is not independent of the standpoint of the observer with the watch or clock, as we know from 
experience”. 

The message by Einstein is simply illustrated in Eig. 6. If we consider in an inertial frame a source (solid line) 
moving with high speed and emitting light signals (dashed lines) along the direction of its motion, a far away observer 
will measure a delay Af^ between the arrival time of two signals emitted at the origin and after a time interval Af in the 
laboratory frame. The real velocity of the source is given by: 


V = 


Ar 

At 


( 1 ) 


and the apparent velocity is given by: 

Ar 

= —, ( 2 ) 

As pointed out by Einstein the adoption of coordinating light signals simply by their arrival time as in Eq.(2), without 
an adequate definition of synchronization, is incorrect and leads to unsurmountable difficulties as well as to apparently 
“superluminal” velocities as soon as motions close to the speed of light are considered. 

The use of Ata as a time coordinate, often tacitly adopted by astronomers, should be done, if at all, with proper care. 
The relation between Ata and the correct time parameterization in the laboratory frame has to be taken into account: 
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In other words, the relation between the arrival time and the laboratory time cannot be done without a knowledge of 
the speed along the entire world-line of the source. In the case of GRBs, such a worldline starts at the moment of 
gravitational collapse. It is of course clear that the parameterization in the laboratory frame has to take into account 
the cosmological redshift z of the source. We then have, at the detector: 

Ati = {\+z)Ata. (4) 


In the current GRB literature, Eq.(3) has been systematically neglected by addressing only the afterglow description 
neglecting the previous history of the source. Often the integral equation has been approximated by a clearly incorrect 
instantaneous value: 

At 

^ • ( 5 ) 

The attitude has been adopted that it should be possible to consider separately the afterglow part of the GRB 
phenomenon, without the knowledge of the entire equation of motion of the source. 

This point of view has reached its most extreme expression in the works reviewed by Piran [29, 30], where the 
so-called “prompt radiation”, lasting on the order of 10^ s, is considered as a burst emitted by the prolonged activity of 
an “inner engine”. In these models, generally referred to as the “internal shock model”, the emission of the afterglow 
is assumed to follow the “prompt radiation” phase [31, 32, 33, 34, 35]. 

As we outline in the following, such an extreme point of view originates from the inability of obtaining the time 
scale of the “prompt radiation” from a burst structure. These authors consequently appeal to the existence of an “ad 
hoc” inner engine in the GRB source to solve this problem. 




Figure 6. Relation between the arrival time ta and the laboratory time t. Details in Ruffini et al. [13, 28]. 


We show in the following how this difficulty has been overcome in our approach by interpreting the “prompt 
radiation” as an integral part of the afterglow and not as a burst. This explanation can be reached only through a 
relativistically correct theoretical description of the entire afterglow (see section ). Within the framework of special 
relativity we show that it is not possible to describe a GRB phenomenon by disregarding the knowledge of the 
entire past worldline of the source. We show that at 10^ seconds the emission occurs from a region of dimensions 
of approximately lO'® cm, well within the region of activity of the afterglow. This point was not appreciated in the 
current literature due to the neglect of the apparent superluminal effects implied by the use of the “pathological” 
parametrization of the GRB phenomenon by the arrival time of light signals. 





The effective potential 
of a Kerr-Newmann geometry 



Figure 7. The effective potential corresponding to the circular orbits in the equatorial plane of a black hole is given as a function 
of the angular momentum of the test particle. This digram was originally derived by Ruffini and Wheeler (right picture). For details 
see Landau and Lifshitz [9] and Rees, Ruffini and Wheeler [40]. 


An additional difference between our treatment and the ones in the current literature relates to the assumption of the 
existence of scaling laws in the afterglow phase: the power law dependence of the Lorentz gamma factor on the radial 
coordinate is usually systematically assumed. From the proper use of the relativistic transformations and by the direct 
numerical and analytic integration of the special relativistic equations of motion we demonstrate (see section ) that 
no simple power-law relation can be derived for the equations of motion of the system. This situation is not new for 
workers in relativistic theories: scaling laws exist in the extreme ultrarelativistic regimes and in the Newtonian ones 
but not in the intermediate fully relativistic regimes (see e.g. Ruffini [36]). 


GRBs and general relativity 

Three of the most important works in the field of general relativity have certainly been the discovery of the Kerr 
solution [37], its generalization to the charged case (Newman et al. [38]) and the formulation by Brandon Carter [39] 
of the Hamilton-Jacobi equations for a charged test particle in the metric and electromagnetic field of a Kerr-Newman 
solution (see e.g. Landau and Lifshitz [9]). The equations of motion, which are generally second order differential 
equations, were reduced by Carter to a set of first order differential equations which were then integrated by using 
an effective potential technique by Ruffini and Wheeler for the Kerr metric (see e.g. Landau and Lifshitz [9]) and by 
Ruffini for the Reissner-Nordstrom geometry (Ruffini [36], see Fig. 7). 




















































All the above mathematical results were essential for understanding the new physics of gravitationally collapsed 
objects and allowed the publication of a very popular article: “Introducing the black hole” (Ruffini and Wheeler [7]). 
In that paper, we advanced the ansatz that the most general black hole is a solution of the Einstein-Maxwell equations, 
asymptotically flat and with a regular horizon: the Kerr-Newman solution, characterized only by three parameters: the 
mass M, the charge Q and the angular momentum L. This ansatz of the “black hole uniqueness theorem” still today 
after thirty years presents challenges to the mathematical aspects of its complete proof (see e.g. Carter [41] and Bini 
et al. [42]). In addition to these mathematical difficulties, in the field of physics this ansatz contains the most profound 
consequences. The fact that, among all the possible highly nonlinear terms characterizing the gravitationally collapsed 
objects, only the ones corresponding solely to the Einstein Maxwell equations survive the formation of the horizon has, 
indeed, extremely profound physical implications. Any departure from such a minimal configuration either collapses 
on the horizon or is radiated away during the collapse process. This ansatz is crucial in identifying precisely the process 
of gravitational collapse leading to the formation of the black hole and the emission of GRBs. Indeed, in this specific 
case, the Born-like nonlinear [43] term of the Heisenberg-Euler-Schwinger [25, 26] Lagrangian are radiated away 
prior to the formation of the horizon of the black hole (see e.g. Ruffini et al. [27]). Only the nonlinearity corresponding 
solely to the classical Einstein-Maxwell theory is left as the outcome of the gravitational collapse process. 

The same effective potential technique (see Landau and Lifshitz [9]), which allowed the analysis of circular orbits 
around the black hole, was crucial in reaching the equally interesting discovery of the reversible and irreversible 
transformations of black holes by Christodoulou and Ruffini [11], which in turn led to the mass-energy formula of the 
black hole: 
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is the horizon surface area, is the irreducible mass, r+ is the horizon radius and p+ is the quasi-spheroidal 
cylindrical coordinate of the horizon evaluated at the equatorial plane. Extreme black holes satisfy the equality in 
Eq.(7). 

Erom Eq.(6) follows that the total energy of the black hole Ebh can be split into three different parts: rest mass. 
Coulomb energy and rotational energy. In principle both Coulomb energy and rotational energy can be extracted from 
the black hole (Christodoulou and Ruffini [11]). The maximum extractable rotational energy is 29% and the maximum 
extractable Coulomb energy is 50% of the total energy, as clearly follows from the upper limit for the existence of a 
black hole, given by Eq.(7). We refer in the following to both these extractable energies as the blackholic energy. 

The existence of the black hole and the basic correctness of the circular orbits has been proven by the observations 
of Cygnus-Xl (see e.g. Giacconi and Ruffini [6]). However, in binary X-ray sources, the black hole uniquely acts 
passively by generating the deep potential well in which the accretion process occurs. It has become tantalizing to look 
for astrophysical objects in order to verify the other fundamental prediction of general relativity that the blackholic 
energy is the largest energy extractable from any physical object. 

As we shall see in the next section, the feasibility of the extraction of the blackholic energy has been made possible 
by the quantum processes of creating, out of classical fields, a plasma of electron-positron pairs in the field of black 
holes. The manifestation of such process of energy extraction from the black hole is astrophysically manifested by the 
occurrence of GRBs. 


GRBs and quantum electro-dynamics 

That a static electromagnetic field stronger than a critical value: 
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can polarize the vacuum and create electron-positron pairs was clearly evidenced by Heisenberg and Euler [25]. The 
major effort in verifying the correctness of this theoretical prediction has been directed in the analysis of heavy ion 
collisions (see Ruffini et al. [27] and references therein). From an order-of-magnitude estimate, it appears that around 
a nucleus with a charge: 

Z,~^~137 (10) 

the electric field can be stronger than the electric field polarizing the vacuum. A more accurate detailed analysis taking 
into account the bound states levels around a nucleus brings to a value of 

Z,~173 (11) 

for the nuclear charge leading to the existence of a critical field. From the Heisenberg uncertainty principle it follows 
that, in order to create a pair, the existence of the critical field should last a time 

At-(12) 

which is much longer then the typical confinement time in heavy ion collisions which is 


At-(13) 

nipC^ 

This is certainly a reason why no evidence for pair creation in heavy ion collisions has been obtained although 
remarkable effort has been spent in various accelerators worldwide. Similar experiments involving laser beams 
encounter analogous difficulties (see e.g. Ruffini et al. [27] and references therein). 

The alternative idea was advanced in 1975 [23] that the critical field condition given in Eq.(9) could be reached 
easily, and for a time much larger than the one given by Eq.(12), in the field of a Kerr-Newman black hole in a range 
of masses 3.2Mq < Mbh < 7.2 x IO^Mq. In that paper we have generalized to the curved Kerr-Newman geometry 
the fundamental theoretical framework developed in Minkowski space by Heisenberg-Euler [25] and Schwinger [26]. 
This result was made possible by the work on the structure of the Kerr-Newman spacetime previously done by Carter 
[39] and by the remarkable mathematical craftsmanship of Thibault Damour then working with me as a post-doc in 
Princeton. 

The maximum energy extractable in such a process of creating a vast amount of electron-positron pairs around a 
black hole is given by: 

E^ax = 1.8 X 10^4 (Mbh/Mq) erg. (14) 

We concluded in that paper that such a process “naturally leads to a most simple model for the explanation of the 
recently discovered y-rays bursts”. 

At that time, GRBs had not yet been optically identified and nothing was known about their distance and conse¬ 
quently about their energetics. Fiterally thousands of theories existed in order to explain them and it was impossible 
to establish a rational dialogue with such an enormous number of alternative theories. We did not pursue further our 
model until the results of the BeppoSAX mission, which clearly pointed to the cosmological origin of GRBs, implying 
for the typical magnitude of their energy precisely the one predicted by our model. 

It is interesting that the idea of using an electron-positron plasma as a basis of a GRB model was independently 
introduced years later in a set of papers by Cavallo and Rees [44], Cavallo and Horstman [45] and Horstman and 
Cavallo [46]. These authors did not address the issue of the physical origin of their energy source. They reach their 
conclusions considering the pair creation and annihilation process occurring in the confinement of a large amount of 
energy in a region of dimension ~ 10 km typical of a neutron star. No relation to the physics of black holes nor to the 
energy extraction process from a black hole was envisaged in their interesting considerations, mainly directed to the 
study of the opacity and the consequent dynamics of such an electron-positron plasma. 

After the discovery of the afterglows and the optical identification of GRBs at cosmological distances, implying 
exactly the energetics predicted in Eq.(14), we returned to the analysis of the vacuum polarization process around a 
black hole and precisely identified the region around the black hole in which the vacuum polarization process and 
the consequent creation of electron-positron pairs occur. We defined this region, using the Greek name dyad for pairs 
(Svag, SvaSog), to be the “dyadosphere” of the black hole, bounded by the black hole horizon and the dyadosphere 
radius given by (see Ruffini [47], Preparata et al. [48] and Fig.8): 
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Figure 8. The dyadosphere is comprised between the horizon radius and the radius of the dyadosphere. All this region is filled 
with electron-positron pairs and photons in thermal equilibrium. Details in Ruffini [47], Preparata et al. [48], Ruffini et al. [49]. 


where we have introduced the dimensionless mass and charge parameters /r = Mbh/Mq, Qj{M bhVG) < 1. 

The analysis of the dyadosphere was developed, at that time, around an already formed black hole. In recent months 
we have been developing the dynamical formation of the black hole and correspondingly of the dyadosphere during 
the process of gravitational collapse, reaching some specific signatures which may be detectable in the structure of the 
short and long GRBs (Cherubini et al. [50], Ruffini and Vitagliano [51, 52], Ruffini et al. [49, 53, 54]). 


THE DYNAMICAL PHASES FOLLOWING THE DYADOSPHERE FORMATION 

Many details of this topic have been presented in great details in Ruffini et al. [28]. 

After the vacuum polarization process around a black hole, one of the topics of the greatest scientific interest is 
the analysis of the dynamics of the electron-positron plasma formed in the dyadosphere. This issue was addressed by 
us in a very effective collaboration with Jim Wilson at Livermore. The numerical simulations of this problem were 
developed at Livermore, while the semi-analytic approach was developed in Rome (Ruffini et al. [55]). 

The corresponding treatment in the framework of the Cavallo et al. analysis was performed by Piran et al. [56] 
also using a numerical approach, by Bisnovaty-Kogan and Murzina [57] using an analytic approach and by Meszaros, 
Laguna and Rees [58] using a numerical and semi-analytic approach. 

Although some analogies exists between these treatments, they are significantly different in the theoretical details 
and in the final results. Since the final result of the GRB model is extremely sensitive to any departure from the correct 
treatment, it is indeed very important to detect at every step the appearance of possible fatal errors. 

A conclusion common to all these treatments is that the electron-positron plasma is initially optically thick and 
expands till transparency reaching very high values of the Lorentz gamma factor. A second point, which is common. 
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Figure 9. The optically thick phase of our model are qualitatively represented in this diagram. There are clearly recognizable 1) 
the PEM pulse phase, 2) the impact on the baryonic remnant, 3) the PEMB pulse phase and the final approach to transparency with 
the emission of the P-GRB (see Eig. 10). Details in Ruffini et al. [28]. 


is the discovery of a new clear feature: the plasma shell expands but the Lorentz contraction is such that its width in 
the laboratory frame appears to be constant. 

There is however a major difference between our approach and the ones of Piran, Meszaros and Rees, in that 
the dyadosphere is assumed by us to be filled uniquely with an electron-positron plasma. Such a plasma expands in 
substantial agreement with the results presented in the work of Bisnovati-Kogan and Murzina [57]. In our model the 
pulse of electron-positron pairs and photons (PEM Pulse, see Ruffini et al. [55]) evolves and at a radius on the order 
of lO'** cm it encounters the remnant of the star progenitor of the newly formed black hole. The PEM pulse is then 
loaded with baryons. A new pulse is formed of electron-positron-photons and baryons (PEMB Pulse, see Ruffini et al. 
[59]) which expands all the way until transparency is reached. At transparency the emitted photons give origin to what 
we define as the Proper-GRB (see Ruffini et al. [14] and Pig. 9). 

In our approach, the baryon loading is measured by a dimensionless quantity 


B = 


Mbc^ 

^dya 


(16) 


which gives direct information about the mass Mb of the remnant. The corresponding treatment done by Piran and 
collaborators (Shemi & Piran [60], Piran et al. [56]) and by Meszaros, Laguna and Rees [58] differs in one important 
respect: the baryonic loading is assumed to occur since the beginning of the electron-positron pair formation and no 
relation to the mass of the remnant of the collapsed progenitor star is attributed to it. 

A marked difference also exists between our description of the rate equation for the electron-positron pairs and 
the ones by those authors. While our results are comparable with the ones obtained by Piran under the same initial 
conditions, the set of approximations adopted by Meszaros, Laguna and Rees [58] appears to be too radical and leads 
to very different results violating energy and momentum conservation (see Bianco et al. [61]). 

Prom our analysis (Ruffini et al. [59]) it also becomes clear that such expanding dynamical evolution can only occur 
for values of Z? < 10^^. This prediction, as we will show shortly in the three GRB sources considered here, is very 
satisfactorily confirmed by observations. 

Prom the value of the B parameter, related to the mass of the remnant, it therefore follows that the collapse to a 
black hole leading to a GRB is drastically different from the collapse to a neutron star. While in the case of a neutron 
star collapse a very large amount of matter is expelled, in many instances well above the mass of the neutron star itself, 
in the case of black holes leading to a GRB only a very small fraction of the initial mass (^ 10^^ or less) is expelled. 
The collapse to a black hole giving rise to a GRB appears to be much smoother than any collapse process considered 
until today: almost 99.9% of the star has to be collapsing simultaneously! 





Detector arrival time (tg'*) (s) 

Figure 10. The P-GRB emitted at the transparency point at a time of arrival which has been computed following the 
prescriptions of Eq.(3). Details in Ruffini et al. [14, 28]. 


We summarize in Figs. 9-10 the optically thick phase of GRBs in our model; we start from a given dyadosphere of 
energy E^ya, the pair-electromagnetic pulse (PEM pulse) self-accelerates outward typically reaching Lorentz gamma 
factors 7 ^ 200 at r ^ 10^** cm; at this point the collision of the PEM pulse with the remnant of the progenitor star 
occurs with an abrupt decrease in the value of the Lorentz gamma factor; a new pair-electromagnetic-baryon pulse 
(PEMB pulse) is formed which self-accelerates outward until the system becomes transparent. 

The photon emission at this transparency point is the Proper-GRB (P-GRB). An accelerated beam of baryons with 
an initial Lorentz gamma factor jo starts to interact with the interstellar medium at typical distances from the black 
hole of To ^ lO*'* cm and at a photon arTival time at the detector on the Earth surface of ~ 0.1 s. These values 
determine the initial conditions of the afterglow. 


THE DESCRIPTION OF THE AFTERGLOW 

After reaching transparency and the emission of the P-GRB, the accelerated baryonic matter (the ABM pulse) interacts 
with the interstellar medium (ISM) and gives rise to the afterglow (see Pig. 11). Also in the descriptions of this last 
phase many differences exist between our treatment and the other ones in the current literature. 
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Figure 11. The GRB afterglow phase is here represented together with the optically thick phase (see Fig. 9). The value of the 
Lorentz gamma factor is here given from the transparency point all the way to the ultrarelativisitc, relativistic and non relativistic 
regimes. Details in Ruffini et al. [28]. 


The initial value problem 

The initial conditions of the afterglow era are determined at the end of the optically thick era when the P-GRB is 
emitted. As recalled in the last section, the transparency condition is determined by a time of arrival a value of 
the gamma Lorentz factor jo, a value of the radial coordinate ro, an amount of baryonic matter Mb which are only 
functions of the two parameters E^iya and B (see Eq.(16)). It is appropriate here to emphasize again that, in order to 
have the expansion leading to an observed GRB, one must have B < 10^^. 

This connection to the optically thick era is missing in the current approach in the literature which attributes the 
origin of the “prompt radiation” to an unspecified inner engine activity (see Piran [29] and references therein). The 
initial conditions at the beginning of the afterglow era are obtained by a best fit of the later parts of the afterglow. This 
approach is quite unsatisfactory since, as we will explicitly show, the theoretical treatments currently adopted in the 
description of the afterglow are not correct. The fit using an incorrect theoretical treatment leads necessarily to the 
wrong conclusions as well as, in turn, to the determination of incorrect initial conditions. 


The equations of the afterglow dynamics 

Let us first summarize the commonalities between our approach and the ones in the current literature. In both cases 
(see Piran [29], Chiang & Dermer [62] and Ruffini et al. [28]) a thin shell approximation is used to describe the 
collision between the ABM pulse and the ISM; 


^E'mt — 

(y , 

(17a) 

dj = 


(17b) 

dM = 

~p~dE[j^i H" dM'i^xn ) 

(17c) 

dMism — 

Anmpnismr^dr, 

(17d) 


where isint, 7 and M are respectively the internal energy, the Lorentz factor and the mass-energy of the expanding 
pulse, tiism is the ISM number density which is assumed to be constant, liip is the proton mass, e is the emitted fraction 
of the energy developed in the collision with the ISM and Mi^m is the amount of ISM mass swept up within the radius 
r: Mism = (4/3);r(r^ — ro^)mpnism, where ro is the starting radius of the shock front. In general, an additional condition 
is needed in order to determine e as a function of the radial coordinate. In the following, e is assumed to be constant 
and such an approximation appears to be correct in the GRB context. 







In both our work and in the current literature (see Piran [29], Chiang & Dermer [62] and Ruffini et al. [28]) a first 
integral of these equations has been found, leading to expressions for the Lorentz gamma factor as a function of the 
radial coordinate. In the “fully adiabatic condition” (i.e. e = 0) we have: 


7^ = 


yl + 2yo (M;,^IMb) + 
1 + 27 o {M;,JMb) + 


while in the “fully radiative condition” (i.e. e = 1 ) we have: 


7 = 


1 + (1 + y-‘) [1 + (1/2) (Mism/Mfi)] 

Jo^ + {Mis JMb) (1 + 770 [1 + (1/2) (M.^/Mb)] ’ 


(18) 


(19) 


where jo and Mb are respectively the values of the Lorentz gamma factor and of the mass of the accelerated baryons 
at the beginning of the afterglow phase and ro is the value of the radius r at the beginning of the afterglow phase. 

A major difference between our treatment and the other ones in the current literature is that we have integrated the 
above equations analytically, obtaining the explicit analytic form of the equations of motion for the expanding shell 
in the afterglow for a constant ISM density. For the fully radiative case we have explicitly integrated the differential 
equation for r (f) in Eq.(19), recalling that 7 ^^ = 1 — [dr/ (cdt)]^, where t is the time in the laboratory reference frame. 
We have then obtained a new explicit analytic solution of the equations of motion for the relativistic shell in the entire 
range from the ultra-relativistic to the non-relativistic regimes: 


. MB-m° J [A+{r/--o)]3(A3 + l) ^ mjr, 

+ -- 
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[A3 + (r/r„)^'](A+l)^' j &cVC 




( 20 ) 
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arctan — arctan 
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where A = /J{Mb - m°) /m°, C = Mb^{Jo - l)/(7o + 1) and m° = (4/3) 7rmp,Hsm'o 
Correspondingly, in the adiabatic case we have: 


■ to ’ 


( 21 ) 


In the current literature, following Blandford and McKee [63], a so-called “ultrarelativistic” approximation jo ^ 
7^1 has been widely adopted by many authors to solve Eqs.(17) (see e.g. Sari [64, 65], Waxman [ 66 ], Rees & 
Meszaros [67], Granot et al. [ 68 ], Panaitescu & Meszaros [69], Piran [29], Gruzinov & Waxman [70], van Paradijs et 
al. [71], Meszaros [72] and references therein). This leads to simple constant-index power-law relations: 


70c r 


( 22 a) 


with a = 3 in the fully radiative case and a = 3/2 in the fully adiabatic case. This simple relation is in stark contrast 
to the complexity of Eq.(19) and Eq.(18) respectively. In the same spirit, instead of Eq.(20) and Eq.(21), some authors 
have assumed the following much simpler approximation for the relation between the time and the radial coordinate 
of the expanding shell, both in the fully radiative and in the fully adiabatic cases: 


ct = r, 

while others, like e.g. Panaitescu & Meszaros [69], have integrated the approximate Eq.(22a), obtaining: 

cf = r [l-f (4fl-|-2)^' 7 ^^ (r)j . 


( 22 b) 


( 22 c) 


Again, it is appropriate here to emphasize the stark contrast between Eqs.(22b),(22c) and the exact analytic solutions 
of Eqs.(17), expressed in Eqs.(20),(21). 


The equitemporal surfaces (EQTSs) 

As pointed out long ago by Couderc [73], in all relativistic expansion the crucial geometrical quantities with respect 
to a physical observer are the “equitemporal surfaces” (EQTSs), namely the locus of source points of the signals 
arriving at the observer at the same time. 




Figure 12. Not all values of iJ are allowed. Only photons emitted at an angle such that cost? > (v/c) can be viewed by the 
observer. Thus the maximum allowed value corresponds to cos = (r/c). In this figure we show j}max (i-e. the angular 
amplitude of the visible area of the ABM pulse) in degrees as a function of the arrival time at the detector for the photons emitted 
along the line of sight (see text). In the earliest GRB phases v ~ c and so i}max ~ 0. On the other hand, in the latest phases of 
the afterglow the ABM pulse velocity decreases and j}max tends to the maximum possible value, i.e. 90°. Details in Ruffini et al. 
[76, 28] 


For a relativistically expanding spherically symmetric source the EQTSs are surfaces of revolution about the line of 
sight. The general expression for their profile, in the form tJ = tJ(r), corresponding to an arrival time G of the photons 
at the detector, can be obtained from (see e.g. Ruffini et al. [28], Bianco and Ruffini [74, 75] and Figs. 12-14): 

eta = ct (r) — rCOS'/}+ r* , (23) 

where r* is the initial size of the expanding source, 1 } is the angle between the radial expansion velocity of a point on 
its surface and the line of sight, and f =t{r) is its equation of motion, expressed in the laboratory frame, obtained by 
the integration of Eqs.(17). Erom the definition of the Lorentz gamma factor 7 ^^ = 1 — (dr/edt)^, we have in fact: 

= ^ [ 1 - 7 ^^ ('■')] (24) 

where 7 (r) comes from the integration of Eqs.(17). 





ABM pulse radius (r) (cm) 

Figure 13. The diameter of the visible area is represented as a function of the ABM pulse radius. In the earliest expansion phases 
(y ~ 310) T^mojc is very small (see left pane and Fig. 14), so the visible area is just a small fraction of the total ABM pulse surface. 
On the other hand, in the final expansion phases 90° and almost all the ABM pulse surface becomes visible. Details in 

Ruffini et al. [76, 28] 






Figure 14. This figure shows the temporal evolution of the visihle area of the ABM pulse. The dashed half-circles are the 
expanding ABM pulse at radii corresponding to different laboratory times. The black curve marks the boundary of the visible 
region. The black hole is located at position (0,0) in this plot. Again, in the earliest GRB phases the visible region is squeezed along 
the line of sight, while in the final part of the afterglow phase almost all the emitted photons reach the observer. This time evolution 
of the visible area is crucial to the explanation of the GRB temporal structure. Details in Ruffini et al. [76, 28] 


The two EQTSs are represented at selected values of the arrival time G in Fig. 16, where the illustrative case of 
GRB 991216 has been used as a prototype. The initial conditions at the beginning of the afterglow era are in this 
case given by jo = 310.131, ro = 1.943 x 10*^ cm, G = 6.481 x 10^ s, r* = 2.354 x 10^ cm (see Ruffini et al. 
[13, 14, 76, 28]). 


The bolometric luminosity of the source 


We assume that the internal energy due to kinetic collision is instantly radiated away and that the corresponding 
emission is isotropic. As in section , let Ae be the internal energy density developed in the collision. In the comoving 
frame the energy per unit of volume and per solid angle is simply 


r dE \ _Ae 
\dVdD.)^ 47t 


(27) 


due to the fact that the emission is isotropic in this frame. The total number of photons emitted is an invariant quantity 
independent of the frame used. Thus we can compute this quantity as seen by an observer in the comoving frame 
(which we denote with the subscript “o”) and by an observer in the laboratory frame (which we denote with no 
subscripts). Doing this we find: 
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dtdQ.d'L 


dtdO.d'L 


A ^cosrJ, 


(28) 
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Figure 15. Due to the extremely high and extremely varying Lorentz gamma factor, photons reaching the detector on the Earth at 
the same arrival time are actually emitted at very different times and positions. We represent here the surfaces of photon emission 
corresponding to selected values of the photon arrival time at the detector: the equitemporal surfaces (EQTS). Such surfaces differ 
from the ellipsoids described by Rees in the context of the expanding radio sources with typical Lorentz factor y ~ 4 and constant. 
In fact, in GRB 991216 the Lorentz gamma factor ranges from 310 to 1. The EQTSs represented here (solid lines) correspond 
respectively to values of the arrival time ranging from 5 j (the smallest surface on the left of the plot) to 60 i (the largest one on the 
right). Each surface differs from the previous one by 5 j. To each EQTS contributes emission processes occurring at different values 
of the Lorentz gamma factor. The dashed lines are the boundaries of the visible area of the ABM pulse and the black hole is located 
at position (0,0) in this plot. Note the different scales on the two axes, indicating the very high EQTS “effective eccentricity”. The 
time interval from 5s to 60s has been chosen to encompass the E-APE emission, ranging from y = 308.8 to y = 56.84. Details in 
Ruffini et al. [76, 28] 


where cos tJ comes from the projection of the elementary surface of the shell on the direction of propagation and 
A = 7(1 — j3 cos tJ) is the Doppler factor introduced in the two following differential transformation 

dO-o =dQ.x A-^ (29) 


for the solid angle transformation and 


dto = dt X A * 


(30) 


for the time transformation. The integration in c/E is performed over the visible area of the ABM pulse at laboratory 
time t, namely with 0 < tJ < Td-max and d-max defined in section (see Figs. 12-14). An extra A factor comes from the 
energy transformation: 


Eo = E X A. 


( 31 ) 






Figure 16. Comparison between EQTSs in the adiabatic regime (solid lines) and in the fully radiative regime (dashed lines). The 
left plot shows the EQTSs for ta = 5 s,ta = 15 s, = 30 s and ta = 45 s, respectively from the inner to the outer one. The right plot 
shows the EQTS at an arrival time of 2 days. Details in Bianco and Ruffini [75]. 


See also Chiang and Dermer [62]. Thus finally we obtain: 

dE f dE \ 

-= - A cos V 

dldEidY. \dtdQ.dT, J ^ 


(32) 


Doing this we clearly identify as the energy density in the comoving frame up to a factor ^ (see Eq.(27)). 

Then we have: 


dE 

dtdEl 


L 


— V cos tJ A ^ c/E, 

shell TTT 


(33) 


where the integration in c/E is performed over the ABM pulse visible area at laboratory time f, namely with 0 < < 

Tdmax and T&max defined in section . Eq.(33) gives us the energy emitted toward the observer per unit solid angle and per 
unit laboratory time t in the laboratory frame. 

What we really need is the energy emitted per unit solid angle and per unit detector arrival time so we must use 
the complete relation between and t given in Eq.(23). Eirst we have to multiply the integrand in Eq.(33) by the factor 
(c/f /dtf) to transform the energy density generated per unit of laboratory time t into the energy density generated per 
unit arrival time t^. Then we have to integrate with respect to c/E over the equitemporal surface (EQTS, see section 
) of constant arrival time instead of the ABM pulse visible area at laboratory time t. The analog of Eq.(33) for the 
source luminosity in detector arrival time is then: 


—— V cos t> A —tc/E 

dt^dO. Jeqts 47r dtfj 


(34) 


It is important to note that, in the present case of GRB 991216, the Doppler factor A^^ in Eq.(34) enhances the apparent 
luminosity of the burst, as compared to the intrinsic luminosity, by a factor which at the peak of the afterglow is in the 
range between lO'® and lO'^! 

We are now able to reproduce in Eig. 17 the general behavior of the luminosity starting from the P-GRB to the 
latest phases of the afterglow as a function of the arrival time. It is generally agreed that the GRB afterglow originates 
from an ultrarelativistic shell of baryons with an initial Lorentz factor /o ~ 200-300 with respect to the interstellar 
medium (see e.g. Ruffini et al. [28], Bianco & Ruffini [74] and references therein). Using GRB 991216 as a prototype, 
in Ruffini et al. [13, 14] we have shown how from the time varying bolometric intensity of the afterglow it is possible 
to infer the average density (n,jm) = 1 particle/cm^ of the Interstellar Medium (ISM) in a region of approximately 
10^^ cm surrounding the black hole giving rise to the GRB phenomenon. 

It was shown in Ruffini et al. [76] that the theoretical interpretation of the intensity variations in the prompt phase 
in the afterglow implies the presence in the ISM of inhomogeneities of typical scale 10*^ cm. Such inhomogeneities 
were there represented for simplicity as spherically symmetric over-dense regions with ~ 10^ {riism) separated 

by under-dense regions with ~ 10^^ (n,im) also of typical scale ~ 10^^ cm in order to keep (n,jm) constant. 










Figure 17. Bolometric luminosity of P-GRB and afterglow as a function of the arrival time. Details in Ruffini et al. [28]. 
Reproduced and adapted from Ruffini et al. [77] with the kind permission of the publisher. 


The summary of these general results are shown in Fig. 18, where the P-GRB, the emission at the peak of the 
afterglow in relation to the “prompt emission” and the latest part of the afterglow are clearly identified for the source 
GRB 991216. Details in Ruffini et al. [28]. 


THE THEORY OF THE LUMINOSITY IN FIXED ENERGY BANDS AND SPECTRA 

OF THE AFTERGLOW 

Having obtained a general agreement between the observed luminosity variability and our treatment of the bolometric 
luminosity, we have further developed the model in order to explain 

a) the details of the observed luminosity in fixed energy bands, which are the ones actually measured by the detectors 
on the satellites 

b) the instantaneous as well as the average spectral distribution in the entire afterglow and 

c) the observed hard to soft drift observed in GRB spectra. 

In order to do so we have developed (Ruffini et al. [78]) a more detailed theory of the structure of the shock front 
giving rise to the afterglow. We have modeled the interaction between the ultrarelativistic shell of baryons and the ISM 
by a shock front with three well-defined layers (see e.g. secs. 85-89, 135 of Landau & Lifshitz [79], ch. 2 and sec. 
13-15 of Zel’dovich & Rayzer [80] and sec. IV, 11-13 of Sedov [81]). From the back end to the leading edge of this 
shock front there is: 
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Figure 18. The detailed features of GRB 991216 evidenced by our theoretical models are here reproduced. The P-GRB, the 
“prompt radiation” and what is generally called the afterglow. It is clear that the prompt radiation coincides with the extended 
afterglow peak emission (E-APE) and has been considered as a burst only as a consequence of the high noise threshold in the 
observations. Details in Ruffini et al. [76, 28]. 


a) A compressed high-temperature layer, of thickness A', in front of the relativistic baryonic shell, created by the 
accumulated material swept up in the ISM. 

b) A thin shock front, with a jump AT in the temperature which has been traditionally estimated in the comoving 
frame by the Rankine-Hugoniot adiabatic equations: 


AT - (3/16)mp5vV/t~ 1.5 x 10“ 


5v/{l0^kms 


K, 


(35) 


where 5v is the velocity jump, irip is the proton mass and k is Boltzmann’s constant. Of course such a treatment, valid 
for 7 ~ 1, has to be modified (see below) in our novel treatment for the 7 ~ 200 case relevant to GRBs. 

c) A pre-shock layer of ISM swept-up matter at much lower density and temperature, both of which change abruptly 
at the thin shock front behind it. 

At larger distances ahead of the expanding fireball the ISM is at still smaller densities. The upper limit to the 
temperature jump at the thin shock front, given in Eq.(35), is due to the transformation of kinetic energy to thermal 
energy, since the particle mean free path is assumed to be less then the thickness of the layer (a). The thermal emission 
of the observed X- and gamma ray radiation, which as seen from the observations reveals a high level of stability, is 
emitted in the above region (a) due to the sharp temperature gradient at the thin shock front described in the above 
region (b). 















The optical and radio emission comes in our model from the extended region (c). The description of such a 
region, unlike the sharp and well-defined temperature gradient occurring in region (b), requires magnetohydrodynamic 
simulations of the evolution of the electron energy distribution of the synchrotron emission. Such analysis has been 
performed using 3-D Eulerian MHD codes for the particle acceleration models to produce the energy spectrum of 
cosmic rays at supernova envelope fronts (see e.g. McKee and Cowie [82], Tenorio-Tagle et al. [83], Stone and 
Norman [84], Jun & Jones [85]). Other challenges are the magnetic field and the instabilities. We mention two key 
phenomena: first, the importance of the development of Kelvin-Helmholtz and Rayleigh-Taylor instabilities ahead of 
the thin shock front. The second is the dual effect that the shock front has on the ISM initial magnetic field, first 
through the compression of the swept-up matter containing the field and secondly the amplification of the radial 
magnetic field component due to the Rayleigh-Taylor instability. Simulations of both effects (see e.g. Jun and Jones 
[85] and references therein), modeling the synchrotron radio emission for an expanding supernova shell at various 
initial magnetic field and ISM parameter values, shows for example that the presence of an initial tangential magnetic 
field component may essentially affect the resulting magnetic field configuration and hence the outgoing radio flux 
and spectrum. Among the additional effects to be taken into account are the initial inhomogeneity of the ISM and the 
contribution of magnetohydrodynamic turbulence. 

In our approach we focus uniquely on the X- and gamma ray radiation, which appears to be conceptually much 
simpler than the optical and radio emission. It is perfectly predictable by a set of constitutive equations (see next 
section), which leads to directly verifiable and very stable features in the spectral distribution of the observed GRB 
afterglows. In line with the observations of GRB 991216 and other GRB sources, we assume in the following that the 
X- and gamma ray luminosity represents approximately 90% of the energy flux of the afterglow, while the optical and 
radio emission represents only the remaining 10%. 

This approach differs significantly from the other ones in the current literature, where attempts are made to explain 
at once all the multi-wavelength emission in the radio, optical, X and gamma ray as coming from a common origin 
which is linked to boosted synchrotron emission. Such an approach has been shown to have a variety of difficulties 
(Ghirlanda et al. [86], Preece et al. [87]) and cannot anyway have the instantaneous variability needed to explain the 
structure in the “prompt radiation” in an external shock scenario, which is indeed confirmed by our model. 


The equations determining the luminosity in fixed energy bands 


Here the fundamental new assumption is adopted (see also Ruffini et al. [88]) that the X- and gamma ray radiation 
during the entire afterglow phase has a thermal spectrum in the co-moving frame. The temperature is then given by: 


Ts = [AEint/ {4nr^Axa^)] , 


(36) 


where AZtint is the internal energy developed in the collision with the ISM in a time interval At in the co-moving frame, 
a is the Stefan-Boltzmann constant and 

^%=AefflA, (37) 

is the ratio between the “effective emitting area” of the afterglow and the surface area of radius r. In GRB 991216 
such a factor is observed to be decreasing during the afterglow between: 3.01 x 10^^ 5.01 X 10^(Ruffini et 

al. [88]). 

The temperature in the comoving frame corresponding to the density distribution described in Ruffini et al. [76] is 
shown in Fig. 19. 

We are now ready to evaluate the source luminosity in a given energy band. The source luminosity at a detector 
arrival time f^, per unit solid angle dQ. and in the energy band [vi, V2] is given by (see Ruffini et al. [28, 88]): 


d_E^ 

dt^dQ. 


Ae „ -A dt 
-—V cost? A -r^W (vi, 


(38) 


where Ae = AEmt jV is the energy density released in the interaction of the ABM pulse with the ISM inhomogeneities 
measured in the comoving frame, A = 7(1 — (v/c) cost?) is the Doppler factor, W (vi, V 2 ,Tarr) is an “effective weight” 
required to evaluate only the contributions in the energy band [vi, V2], dT. is the surface element of the EQTS at detector 
arrival time on which the integration is performed (see also Ruffini et al. [76]) and Tarr is the observed temperature 
of the radiation emitted from dH: 


d'arr — Ts/ [7(1 - (v/c)costJ) (1 -l-z)] . 


(39) 




r (cm) 

Figure 19. The temperature in the comoving frame of the shock front corresponding to the density distribution with the six spikes 
A,B,C,D,E,F presented in Ruffini et al.^. The dashed line corresponds to an homogeneous distribution with = 1. Details in 
Ruffini et al. [78]. 


The “effective weight” W (vi, V 2 ,Tarr) is given by the ratio of the integral over the given energy band of a Planckian 
distribution at a temperature Tarr to the total integral aT^^/. 


1 /■''2 ^hvV 

W {Vi,V2,Tarr) = P {Tarr,V) d [ — ] 

aT^rr Jvi \C J 


where p {Tarr, v) is the Planckian distribution at temperature Tarr- 


{Tarr, V) = {2/h^) hV/ - l) 


ON THE TIME INTEGRATED SPECTRA AND THE HARD-TO-SOFT SPECTRAL 

TRANSITION 

We turn now to the much debated issue of the origin of the observed hard-to-soft spectral transition during the 
GRB observations (see e.g. Frontera et al. [89], Ghirlanda et al. [86], Piran [29], Piro et al. [90]). We consider the 
instantaneous spectral distribution of the observed radiation for three different EQTSs: 

• = 10 s, in the early radiation phase near the peak of the luminosity, 

• = 1-45 X 10^ s, in the last observation of the afterglow by the Chandra satellite, and 





Figure 20. The instantaneous spectra of the radiation observed in GRB 991216 at three different EQTS respectively, from top 
to bottom, for = 10 s, = 10“* s and = 1.45 x 10^ s. These diagrams have been computed assuming a constant {tiism) — 1 
particle/cm^ and clearly explains the often quoted hard-to-soft spectral evolution in GRBs. Details in Ruflini et al. [88]. 


• = 10“^ s, chosen in between the other two (see Fig. 20). 

The observed hard-to-soft spectral transition is then explained and traced back to: 

1 . a time decreasing temperature of the thermal spectrum measured in the comoving frame, 

2. the GRB equations of motion, 

3. the corresponding infinite set of relativistic transformations. 

A clear signature of our model is the existence of a common low-energy behavior of the instantaneous spectrum 
represented by a power-law with index a — -1-0.9. This prediction will be possibly verified in future observations. 

Starting from these instantaneous values, we integrate the spectra in arrival time obtaining what is usually fit in 
the literature by the “Band relation” (Band et al. [91]). Indeed we find for our integrated spectra a low energy spectral 
index a = — 1.05 and an high energy spectral index j3 < —16 when interpreted within the framework of a Band relation 
(see Fig. 21). This theoretical result can be submitted to a direct confrontation with the observations of GRB 991216 
and, most importantly, the entire theoretical framework which we have developed can now be applied to any GRB 
source. The theoretical predictions on the luminosity in fixed energy bands so obtained can be then straightforwardly 
confronted with the observational data. 


THE THREE PARADIGMS FOR THE INTERPRETATION OF GRBS 

Having outlined the main features of our model and shown its application to GRB 991216 used as a prototype, before 
addressing the two new sources which are going to be the focus of this presentation, we recall the three paradigms for 
the interpretation of GRBs we had previously introduced. 

The first paradigm, the relative space-time transformation (RSTT) paradigm (Ruffini et al. [13]) emphasizes the 
importance of a global analysis of the GRB phenomenon encompassing both the optically thick and the afterglow 
phases. Since all the data are received in the detector arrival time it is essential to know the equations of motion of all 
relativistic phases with 7 > 1 of the GRB sources in order to reconstruct the time coordinate in the laboratory frame. 





Figure 21. The time-integrated spectrum of the radiation observed in GRB 991216. The low energy part of the curve below 10 
keV is ht by a power-law with index a = —1.05 and the high energy part above 500 keV is ht by a power-law with an index 
< —16. Details in Ruffini et al. [88]. 


see Eq.(3). Contrary to other phenomena in nonrelativistic physics or astrophysics, where every phase can be examined 
separately from the others, in the case of GRBs all the phases are inter-related by their signals received in arrival time 
t^. There is the need, in order to describe the physics of the source, to derive the laboratory time f as a function of the 
arrival time along the entire past worldline of the source using Eq.(4). 

The second paradigm, the interpretation of the burst structure (IBS) paradigm (Ruffini et al. [14]) covers three 
fundamental issues: 

a) the existence, in the general GRB, of two different components: the P-GRB and the afterglow related by precise 
equations determining their relative amplitude and temporal sequence (see Ruffini et al. [28]); 

b) what in the literature has been addressed as the “prompt emission” and considered as a burst, in our model is not a 
burst at all — instead it is just the emission from the peak of the afterglow (see Eig. 18); 

c) the crucial role of the parameter B in determining the relative amplitude of the P-GRB to the afterglow and 
discriminating between the short and the long bursts (see Pig. 22). Both short and long bursts arise from the same 
physical phenomena: the dyadosphere. The absence of baryonic matter in the remnant leads to the short bursts and no 
afterglow. The presence of baryonic matter with B < 10^^ leads to the afterglow and consequently to its peak emission 
which gives origin to the so-called long bursts. 

The third paradigm, the GRB-Supernova Time Sequence (GSTS) paradigm (Ruffini et al. [15]), deals with the 
relation of the GRB and the associated supernova process, and acquires a special meaning in relation to the sources 
GRB 980425 and GRB 030329 as we will show in the following. 

We now shortly illustrate some consequences of these three paradigms. 


Long bursts are E-APEs 

The order of magnitude estimate usually quoted for the the characteristic time scale to be expected for a burst 
emitted by a GRB at the moment of transparency at the end of the expansion of the optically thick phase is given by 
T ~ GM/c^, which for a IOMq black hole will give ~ 10^^ s. There are reasons today not to take seriously such an 
order of magnitude estimate (see e.g. Ruffini et al. [54]). In any case this time is much shorter then the ones typically 





Figure 22. The energy radiated in the P-GRB (the solid line) and in the afterglow (the dashed line), in units of the total energy of 
the dyadosphere (Fjya), are plotted as functions of the B parameter. 


observed in “prompt radiation” of the long bursts, from a few seconds all the way to 10^ s. In the current literature (see 
e.g. Piran [29] and references therein), in order to explain the “prompt radiation” and overcome the above difficulty 
it has been generally assumed that its origin should be related to a prolonged “inner engine” activity preceding the 
afterglow which is not well identified. 

To us this explanation has always appeared logically inconsistent since there remain to be explained not one but 
two very different mechanisms, independent of each other, of similar and extremely large energetics. This approach 
has generated an additional very negative result: it has distracted everybody working in the field from the earlier very 
interesting work on the optically thick phase of GRBs. 

The way out of this dichotomy in our model is drastically different: 1) indeed the optically thick phase exists, is 
crucial to the GRB phenomenon and terminates with a burst: the P-GRB; 2) the “prompt radiation” follows the P-GRB; 
3) the “prompt radiation” is not a burst: it is actually the temporally extended peak emission of the afterglow (E-APE). 
The observed structures of the prompt radiation can all be traced back to inhomogeneities in the interstellar medium 
(see Eig. 18 and Ruffini et al. [76]). 


Short bursts are P-GRBs 

The fundamental diagram determining the relative intensity of the P-GRB and the afterglow as a function of the 
dimensionless parameter B has been shown in Pig. 22. The underlying machine generating the short and the long 
GRBs is identical: in both cases is the dyadosphere. The main difference relates to the amount of baryonic matter 
engulfed by the electron-positron plasma in their optically thick phase prior to transparency. In the limit of small 




B < 10^^ the intensity of the P-GRB is larger and dominates the afterglow. This corresponds to the short bursts. For 
10^^ < B < 10^^ the afterglow dominates the GRBs and we have the so-called “long bursts”. For B > 10^^ we may 
observe a third class of “bursts”, eventually related to a turbulent process occurring prior to transparency (Ruffini et 
al. [59]). This third family should be characterized by smaller values of the Lorentz gamma factors than in the case of 
the short or long bursts. 


The trigger of multiple gravitational collapses 

The relation between the GRBs and the supernovae is one of the most complex aspects to be addressed by our model, 
which needs the understanding of new fields of general relativistic physics in relation to yet unexplored many-body 
solutions in a substantially new astrophysical scenario. 

As we will show in the two systems GRB980425/SN1998bw and GRB030329/SN2003dh which we are going to 
discuss next, there is in each one the possibility of an astrophysical “triptych”"^ formed by: 

1) the formation of the black hole and the emission of the GRB, 

2 ) the gravitational collapse of an evolved companion star, leading to a supernova, 

3) a clearly identified URCA source whose nature appears to be of the greatest interest. 

This new astrophysical scenario presents new challenges: 

a) The identification of the physical reasons of the instability leading to the gravitational collapse of a ~ IOMq star, 
giving origin to the black hole. Such an implosion must occur radially with negligible mass of the remnant (B < 10^^). 

b) The identification of the physical reasons for the instability leading to the gravitational collapse of an evolved 
companion star, giving origin to the supernova. 

c) The theoretical issues related to the URCA sources, which range today in many possible directions: from the physics 
of black holes, to the physical processes occurring in the expanding supernova remnants, and finally to the very exciting 
possibility that we are observing for the first time a newly born neutron star. The main effort in the next sections is to 
show that the detailed understanding we have reached for the GRB phenomenon and its afterglow allows us to state, 
convincingly, that the URCA source, contrary to what established in the current literature, is not part of the GRB nor 
of its afterglow. 

We will draw in the conclusions some considerations on the possible nature of the URCA sources. 


APPLICATIONS 

We illustrate the application of our GRB model to two different systems, which are quite different in the energetics 
but are both related to supernovae: GRB 980425 and GRB 030329. We will let the gradual theoretical understanding 
of the system to unveil the underlying astrophysical scenario. 


GRB 980425 / SN 1998bw 

Approaches in the current literature have always attempted to explain both the supernova and the GRB as two 
aspects of a single phenomenon assuming that the GRB takes his origin from a specially strong and yet unobserved 
supernova process: a hypernova (see Paczyriski [92], Kulkarni [93], Iwamoto [94]). 

We have taken a very different approach, following Cicero’s classic aphorism “divide et impera”, which was 
adopted as the motto of the Roman empire: “divide and conquer”. In this specific case of GRBs, which are indeed 
a very complex system, we plan to divide and identify the truly independent physical constituents and conquer the 
understanding of the underlying astrophysical process. As we will see, this approach will lead to an unexpected and 
much richer scenario. 

In addition to the source GRB 980425 and the supernova SN1998bw, two X-ray sources have been found by 
BeppoSAX in the error box for the location of GRB 980425: a source SI and a source S2 (Pian et al. [95]), which 


A picture or carving in three panels side by side; esp'. an altarpiece with a central panel and two flanking panels half its size that fold over it 
[Webster’s New collegiate dictionary, G. & C. Merriam Co. (Springfield, Massachussets, U.S.A., 1977)] 



GRB 980425 - SN1998bw 



Figure 23. The “divide et impera” concept applied to the system GRB 980425 / SN 1998bw. Four different components are 
identified: the GRB 980425, the SN 1998bw, and the two sources SI and S2. 


have been traditionally interpreted either as a background source or as a part of the GRB afterglow. See Fig. 23. 
Our approach has been: to first comprehend the entire afterglow of GRB 980425 within our theory. This allows the 
computation of the luminosity in given energy bands, the spectra, the Lorentz gamma factors, and more generally of all 
the dynamical aspects of the source. Having characterized the features of GRB 980425, we can gradually approach the 
remaining part of the scenario, disentangling the GRB observations from those of the supernova and then disentangling 
both the GRB and the supernova observations from those of the sources SI and S2. This leads to a natural identification 
of distinct events and to their autonomous astrophysical characterization. 

Our best fit for GRB 980425 corresponds to E^ya = 1.1 x 10"^^ ergs, B = 1 x 10^^ and the ISM average density 
is found to be (riism) — 0.02 particle/cm^. The plasma temperature and the total number of pairs in the dyadosphere 
are respectively T = 1.028MeV and Ng± = 5.3274 x 10^^. The light curve of the GRB is shown in Figs. 24-25. The 
P-GRB is under the threshold and in the case of this source is not observable (see Ruffini et al. [16], Fraschetti et al. 
[17]). 

The characteristic parameter 3^, defining the filamentary structure of the ISM, monotonically decreases from 
4.81 X lO^*** to 2.65 X 10^'^). The results are given in Fig. 26 where the bolometric luminosity is represented together 
with the optical data of SN1998bw, the source SI and the source S2. It is then clear that GRB 980425 is separated both 
from the supernova data and from the sources SI and S2. 

While the occurrence of the supernova in relation to the GRB has already been discussed within the GRB-Supernova 
Time Sequence (GSTS) paradigm (Ruffini et al. [15]), we like to address here a different fundamental issue: the nature 
of the source SI which we have named, in celebration of the work of Gamow and Shoenberg, URCA-1. It is clear, 









GRB 980425 - SN1998bw 



Figure 24. The bolometric luminosity as a function of the arrival time. The peak of the P-GRB is just below the observational 
noise level. 


from the theoretical predictions of the afterglow luminosity, that the URCA-1 cannot be part of the afterglow (see Figs. 
24, 26). There are three different possibilities for the explanation of such source: 

1) Its possible relation to the black hole formed during the process of gravitational collapse leading to the GRB 
emission. 

2) Its possible relation to emission originating in the early phases of the expansion of the supernova remnant. 

3) The very exciting possibility that for the first time we are observing a newly born neutron star out of the supernova 
phenomenon. 

While some general considerations will be discussed in the conclusions, we would like to stress here the paramount 
importance of following the further time history of URCA-1 and of the source S2. If, as we propose, S2 is a background 
source, its flux should be practically constant in time and this source has nothing to do with the GRB 980425 / 
SN1998bw system. The drastic behavior of the URCA-1 luminosity reported in the talk by Elena Plan in this meeting, 
showing the latest URCA-1 observations by the XMM and Chandra satellites, is crucial for the understanding of the 
nature of this source. Some very qualitative luminosity curves are sketched in Fig. 26, illustrating the possible time 
evolution of URCA-1. They are still very undetermined today due to a lack of attention to these observational data and, 
consequently, to the lack of a detailed theoretical model of the phenomenon. We therefore propose to have a dedicated 
attention to the astrophysical “triptych” GRB 980425 / SN 1998bw / URCA-1. 










The prompt emission of GRB 980425 
and the fit in the EMBH model 



Figure 25. The observation by BeppoSAX of the peak of the afterglow in the 40-700 keV energy band is fitted by our model. 


GRB 030329/SN 2003dh 

We have adopted for our modeling of GRB030329 a spherically symmetric distribution for the source and, as initial 
conditions at f = 10^^* s, an e+-e^-photon neutral plasma lying between the radii ri = 2.9 x 10® cm and r 2 = 9.0 x 10^ 
cm. The temperature of such a plasma is 2.1 MeV, the total energy Etot = 2.1 x 10®^ erg and the total number of pairs 
Ng+g- = 1.1 x 10®^. The baryonic matter component Mg is the second free parameter of our theory; B — 4.8 x 10^^. 
At the emission of the P-GRB, the Lorentz gamma factor is Jo = 183.6 and the radial coordinate is Vo = 5.3 x 10^^ 
cm. The ISM average density is best fit by < >= 1 particle/cm^. The third free parameter of our theory is given 

by 1.1 X 10-^<^<5.0x 10-“. 

We then obtain (see also Bernardini et al. [97]) for the GRB 030329 the luminosities in given energy bands, 
computed in the range 2^00 keV with very high accuracy. Figs. 27-30 shows the results for the luminosities in 
the 30-400 keV and 2-10 keV bands, including the “prompt radiation”. Subsequently, the theoretically predicted 
GRB spectra have been evaluated at selected values of the arrival time (Ruffini et al. [98]). 

The splendid news received the evening before the presentation of this talk is graphically represented by the XMM 
observations shown in Fig. 30. Again, the XMM observations, like the corresponding ones of GRB 980425, occur 
after the decaying part of the afterglow and, in analogy to the one occurring in the system GRB 980425 / SN 1998bw / 
URCA-1, we call this source URCA-2. Further observations by XMM are highly recommended to follow the URCA-2 
temporal evolution. Also in this system we are dealing with an astrophysical “triptych”: GRB 030329 / SN 2003dh / 
URCA-2. 
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Figure 26. The bolometric light curves is reported as well as the BeppoSAX MECS observations in the 2-10 keV band of SI 
and S2 (Plan et al. [95]) and the optical data of SN 1998bw (Iwamoto [96]). Here we also report some very qualitative curves to be 
expected for the URCA-1 luminosity. Details in Ruffini et al. [16]. 


ON THE SHORT GRBS 


By the analysis of the hrst and second BATSE catalogs^ Tavani in 1998 [21] (see Fig. 4) conhrmed the previous results 
by Kouveliotou et al. [20] on the existence of two families of GRBs; the so-called “long-bursts” with a soft spectrum 
and duration Af > 2.5sec and the “short-bursts” with harder spectrum and duration Af < 2.5sec. In 2001 we have 
proposed the theory [14] that both short-bursts and long-bursts originate from the same underlying physical process: 
the vacuum polarization of electromagnetic overcritical gravitational collapse leading to the creation of pairs 

at the expenses of the “blackholic” energy [11]. The difference between the short-bursts and long-bursts in this theory 
is mainly due to the amount of baryonic matter, described by the dimensionless parameter B previously mentioned, 
encountered by the pairs in their relativistic expansion (see Fig. 31). Short-bursts occur in a range of B: 

0<B<10^^ (42) 


and the long-bursts occur for: 


10^^ <B < 10^^ 


Compare and contrast Fig. 11 with Fig. 31. 


(43) 


see http://cossc.gsfc.nasa.gov/batse/ 
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Figure 27. The luminosity in the 2-10 keV and in the 30-400 keV energy hands predicted hy our model are fitted to the data of 
R-XTE (GCN Circ. 1996 [99]) and HETE-2 (GCN Circ. 1997 [100]) respectively. The SN 2003dh optical luminosity is given hy 
the crosses (Hjorth et al. [101]). Details in Bemardini et al. [97]. 


An indirect support of our theory was given by Schmidt [103] who has shown that short-bursts and long-bursts have 
the same isotropic-equivalent characteristic peak luminosity. 

In recent work we have systematically developed the theoretical background of the process of gravitational collapse 
of matter involving an electromagnetic field with field strength higher than the critical value for pair creation 
[50, 51, 52, 53, 104, 49]. The goal has been to clarify the physical nature of the process of extracting the “blackholic” 
energy by the creation of matter pairs [11] and to analyze the electromagnetic radiation emission process during 
the transient dynamical phases of the gravitational collapse, leading to the final formation of the black hole. 

All the considerations presented in the description of the long GRBs were based on a dyadosphere of an already 
formed black hole, presented in section . This approximate treatment is very satisfactory in estimating the general 
dependence of the energy of the P-GRB, the kinetic energy of the baryonic matter pulse generating the afterglow and 
consequently the intensity of the afterglow itself. In particular it is possible to obtain the overall time structure of the 
GRB and especially the time of the release of the P-GRB in respect to the moment of gravitational collapse and its 
relative intensity with respect to the afterglow. If, however, we address the issue of the detailed temporal structure of 
the P-GRB and its detailed spectral distribution, the dynamical considerations on the dyadosphere formation, which 
we are going to present in the following sections, are needed (see also [53]). In turn, this detailed analysis is needed 
in order to describe the general relativistic effects close to the horizon formation. As expressed already in section , 
all general relativistic quantum field theory effects are encoded in the fine structure of the P-GRB. As emphasized in 
section , the only way to differentiate between solutions with same but different black hole mass and charge is to 
observe the P-GRBs in the limit B -^0, namely, to observe the short GRBs (see Fig. 31). 














Figure 28. The details of the theoretical fit of the prompt radiation of GRB 030329 have been reproduced by the filamentary 
structure in the ISM in our model. Details in Bernardini et al. [97], 


SOME PROPAEDEUTIC ANALYSIS FOR THE DYNAMICAL FORMATION OF THE 

BLACK HOLE 

While the formation in time of the dyadosphere is the fundamental phenomena we are interested in, we can get an 
insight on the issue of gravitational collapse of an electrically charged star core studying in details a simplified model, 
namely a thin shell of charged dust. 


On the collapsing charged shell in general relativity 


In [105, 106] it is shown that the problem of a collapsing charged shell in general relativity can be reduced to a set 
of ordinary differential equations. We reconsider here the following relativistic system; a spherical shell of electrically 
charged dust which is moving radially in the Reissner-Nordstrom background of an already formed nonrotating black 
hole of mass M\ and charge Qi, with Q\ < Mi. 

The world surface spanned by the shell divides the space-time into two regions: an internal one and an external 
one . The line element in Schwarzschild like coordinate is [50] 


f —f+dt^+f^^dr^ + r^dO? in 

—f^dt^+fl'^dr^ + r^dO? in 


(44) 


where /+ = 1 — ^ -f /_ = 1 — ^ and f_ and f+ are the Schwarzschild-like time coordinates in and 

respectively. M is the total mass-energy of the system formed by the shell and the black hole, measured by an 
observer at rest at infinity and Q — Qo + Qi is the total charge: sum of the charge Qo of the shell and the charge Qi of 
the internal black hole. 


































Figure 29. The perfect fit of the late part of the afterglow of our theoretical model for the 2-10 keV energy bands. The data refers 
to the R-XTE observations (GCN Circ. 1996 [99]). Details in Bernardini et al. [97]. 


Indicating by R the radius of the shell and by T± its time coordinate, the equations of motion of the shell become 

[51] 



dT± 

dx 



(45) 

(46) 


where Mq is the rest mass of the shell and t is its proper time. Eqs.(45,46) (together with Eq.(44)) completely describe 
a 5-parameter (M, Q, Mq) family of solutions of the Einstein-Maxwell equations. Note that Eqs.(45,46) imply 


that 


>0 


(47) 


holds forR >M+ if Q<M and for/? > Mi + iJm\ -Q\ifQ> M. 

Eor astrophysical applications [53] the trajectory of the shell R — R{T+) is obtained as a function of the time 
coordinate T+ relative to the space-time region In the following we drop the -T index from T+. Erom Eqs.(45,46) 
we have _ 


M _ dR^ _ 
dT dz dT 




( 48 ) 

















GRB 030329 - SN2003dh - URCA-2 



Figure 30. The dotted line represents our theoretically predicted GRB030329 light curve in y-rays (30-400 keV) with the horizon¬ 
tal bar corresponding to the mean peak flux from HETE-2 (GCN Circ. 1997 [100]). The solid line represents the corresponding one 
in X-rays (2-10 keV) with the experimental data obtained by R-XTE (GCN Circ. 1996 [99]). The remaining points refer respectively 
to the optical VET data (Hjorth et al. [101]) of SN2003bw and to the X-ray XMM data (Tiengo et al. [102]) of URCA-2. The dash- 
dotted lines corresponds to qualitative luminosity curves expected for URCA-2. It is interesting to compare and contrast these 
results with the ones for GRB980425/SN1998bw (see Fig. 3 in Ruffini et al. [16]). Details in Ruffini et al. [98]. 


where 


F = u{R)^\-m + ^ 

o = r — 

— 2MaR ’ 

Y — M-Mi 
^ - Mo ■ 


(49) 

(50) 

(51) 


Since we are interested in an imploding shell, only the minus sign case in (48) will be studied. We can give the 
following physical interpretation of F. If M — Mi > Mq, F coincides with the Lorentz y factor of the imploding shell 
at infinity; from Eq.(48) it satisfies 



(52) 


When M — Ml < Mq then there is a turning point R*, defined by ^ 
“effective potential” at R* : 


r=x//-(«*)+Mo-' 


2R* 2R 


Ql 


L = 0. In this case F coincides with the 

I K=K' 

< 1 . 


(53) 











The case of short GRBs 



Figure 31. In our theory the short GRBs originate by the same basic process leading also to the long GRBs (see Fig. 11) in the 
limiting case of no baryonic matter present (S —> 0). The optically thick electron-positron pairs expand to ultrarelativistic Lorentz 
gamma factors reaching, then, transparency. The emission at the transparency point originates the short GRBs. 


The solution of the differential equation (48) is given by; 

The functional form of the integral (54) crucially depends on the degree of the polynomial P [R) = (Q? — F) , which 

is generically two, but in special cases has lower values. We therefore distinguish the following cases; 

1. M = Mq + M] ; = Ml , Q = M\ P (R) is equal to 0, we simply have 


/?(r)= const. (55) 

2. M — Mq +Mi; _ ^2 _ q 2 . Q^f^- p ^ constant, we have 

r = const H-— =\(R + 2M)R 
























3. M = Mq + Ml', — P{R) is a first order polynomial and 
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where i2± = Q. (r±). 

4. M ^ Mo + Ml; P (/?) is a second order polynomial and 
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(58) 


Of particular interest is the time varying electric field Sr = % on the external surface of the shell. In order to 


study the variability of Sr with time it is useful to consider in the tridimensional space of parameters {R,T,Sr) the 
parametric curve ‘S : (ji — X, T — T{X), Sr = . In astrophysical applications [53] we are specially interested 


in the family of solutions such that ^ is 0 when R — o° which implies that F = 1. In Fig. 32 we plot the collapse 
curves in the plane {T,R) for different values of the parameter ^ = ^,0 < ^ < 1. The initial data (7o,Po) are chosen 
so that the integration constant in equation (57) is equal to 0. In all the cases we can follow the details of the approach 
to the horizon which is reached in an infinite Schwarzschild time coordinate. In Fig. 32 we plot the parametric curves 
'S in the space (R,T,Sr) for different values of <^. Again we can follow the exact asymptotic behavior of the curves 
"S, Sr reaching the asymptotic value The detailed knowledge of this asymptotic behavior is of great relevance for 

the observational properties of the black hole formation (see e.g. [51]). 

In the case of a shell falling in a flat background (Mi = 2i = 0) Eq.(45) reduces to 
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(59) 


Introducing the total radial momentum P = Mqu'' = Mq^ of the shell, we can express the kinetic energy of the shell 
as measured by static observers in as T = —Mou^^^—Mo = J P^ +Mq —Mq. Then from equation (59) we have 


M = - 


M; 


^ + ^ + \Ip^+mI=Mo + T-^ + ^ 


2R ^ 2R 


,2 
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(60) 


where we choose the positive root solution due to the constraint (47). Eq.(60) is the mass formula of the shell, which 
depends on the time-dependent radial coordinate R and kinetic energy F. If M > Q, a black hole is formed and we 
have 

m 2 , g 2 


M = Mo + T+-^ + §^ 


(61) 


where Tp = T (r+) and = M + a/M^ — is the radius of external horizon of the black hole. 



(S) i 




Figure 32. Left) Collapse curves in the plane (T, R) for M = 20Mq and for different values of the parameter ^. The asymptotic 
behavior is the clear manifestation of general relativistic effects as the horizon of the hlack hole is approached. Right) Electric field 
behavior at the surface of the shell for M = 20Mq and for different values of the parameter ^. The asymptotic behavior is the clear 
manifestation of general relativistic effects as the horizon of the black hole is approached. 


On the physical origin of the terms in mass formula of the hlack hole 

We know from the Christodoulou-Ruffini black hole mass formula that 

M = + (62) 

so it follows that 

Mi„=Mo-^ + r+, (63) 

namely that M-i„ is the sum of only three contributions: the rest mass Mq, the gravitational potential energy and the 
kinetic energy of the rest mass evaluated at the horizon. Mi„ is independent of the electromagnetic energy, a fact noticed 
by Bekenstein [107]. We have taken one further step here by identifying the independent physical contributions XoMi„. 




Next we consider the physical interpretation of the electromagnetic term which can be obtained by evaluating 
the conserved Killing integral 


= j r^dr d cose 
- e! 

“ 2R ’ 


(64) 


where 1.^ is the space-like hypersurface in described by the equation t+ =t = const, with dT/ as its surface ele¬ 
ment vector and where ~ Ppo) the energy-momentum tensor of the electromagnetic 

field. The quantity in Eq.(64) differs from the purely electromagnetic energy 



I rp (em) jy V 



where = f_^_ is the unit normal to the integration hypersurface and grr = /+• This is similar to the analogous 
situation for the total energy of a static spherical star of energy density e within a radius R, m (R) — 4?: Jq dr r^e, which 
differs from the pure matter energy nip (R) = 4n Jq dr^/^r^e by the gravitational energy (see [108]). Therefore the 
term ^ in the mass formula (60) is the total energy of the electromagnetic field and includes its own gravitational 
binding energy. This energy is stored throughout the region 1.^, extending from R to infinity. 


On the energy extraction process of blackholic energy 


We now turn to the problem of extracting the blackholic energy from a black hole (see [11]). We can distinguish 
between two conceptually physically different processes, depending on whether the electric field strength ^ is 

2 3 

smaller or greater than the critical value Here nig and e are the mass and the charge of the electron. As 

already mentioned in this paper an electric field S’ > £[; polarizes the vacuum creating electron-positron pairs (see 

[25]). The maximum value of the electric field around a black hole is reached at the horizon. We then have 

'+ 

the following: 


1. For S+ < Sc the leading energy extraction mechanism consists of a sequence of discrete elementary decay 
processes of a particle into two oppositely charged particles. The condition S-^ < Sc implies 
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(65) 


where Ac is the Compton wavelength of the electron. [109] and [110] have defined as the effective ergosphere the 
region around a black hole where the energy extraction processes occur. This region extends from the horizon r+ 
up to a radius 


^Eerg — 
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( 66 ) 


The energy extraction occurs in a finite number Apo of such discrete elementary processes, each one correspond¬ 
ing to a decrease of the black hole charge. We have 
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(67) 


Since the total extracted energy is (see Eq.(62)) 


(^)pD ~ 
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we obtain for the mean energy per accelerated particle 
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which gives 
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(69) 


One of the crucial aspects of the energy extraction process from a black hole is its back reaction on the irreducible 
mass expressed in [11]. Although the energy extraction processes can occur in the entire effective ergosphere 
defined by Eq. (66), only the limiting processes occurring on the horizon with zero kinetic energy can reach the 
maximum efficiency while approaching the condition of total reversibility (see Fig. 2 in [11] for details). The 
farther from the horizon that a decay occurs, the more it increases the irreducible mass and loses efficiency. Only 
in the complete reversibility limit [11] can the energy extraction process from an extreme black hole reach the 
upper value of 50% of the total black hole energy. 

2. For the leading extraction process is a collective process based on an electron-positron plasma generated 

by the vacuum polarization, (see Fig. 8) as discussed in section III in [28] The condition <#’+ > implies 


GM/c^ 





(70) 


This vacuum polarization process can occur only for a black hole with mass smaller than 2 • 10®M©. The electron- 
positron pairs are now produced in the dyadosphere of the black hole, (note that the dyadosphere is a subregion 
of the effective ergosphere) whose radius is given in Eq.(15). We have r^s 'C reerg- The number of particles 
created and the total energy stored in dyadosphere are given in Eqs.(17,18) of Ref. [51] respectively and we have 
approximately 


(71) 

(72) 

Ea 

The mean energy per particle produced in the dyadosphere {E)^^ = is then 

e~ 

(73) 

which can be also rewritten as 

(£)ds h (^) - ^^IQ^keV. (74) 

Such a process of vacuum polarization, occurring not at the horizon but in the extended dyadosphere region 
{r+ < r < rds) around a black hole, has been observed to reach the maximum efficiency limit of 50% of the total 
mass-energy of an extreme black hole (see e.g. [48]). The conceptual justification of this result follows from the 
present work: the creation process occurs at the expense of the Coulomb energy given by Eq. (64) and 
does not affect the irreducible mass given by Eq. (63), which indeed, as we have proved, does not depend of 
the electromagnetic energy. In this sense, dMi„ = 0 and the transformation is fully reversible. This result will be 
further validated by the study of the dynamical formation of the dyadosphere, which we have obtained using the 
present work and [50] (see [49, 53]). 


Let us now compare and contrast these two processes. 

We have 
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Moreover we see (Eqs. (69), (74)) that (£’)pj 3 is in the range of energies of UHECR, while for ^ ~ 0.1 and M ~ lOM©, 
(Fjjjj, is in the gamma ray range. In other words, the discrete particle decay process involves a small number of particles 
with ultra high energies (~ lO^^eV), while vacuum polarization involves a much larger number of particles with lower 
mean energies (~ IQMeV). 




Figure 33. Space-time diagram of the collapse process leading to the formation of the dyadosphere. As the collapsing core crosses 
the dyadosphere radius the pair creation process starts, and the pairs thermalize in a neutral plasma configuration. Then also the 
horizon is crossed and the singularity is formed. 


Having so established and clarified the basic conceptual processes of the energetic of the black hole, we are now 
ready to approach, using the new analytic solution obtained, the dynamical process of vacuum polarization occurring 
during the formation of a black hole as qualitatively represented in Fig. 33. The study of the dyadosphere dynamical 
formation as well as of the electron-positron plasma dynamical evolution will lead to the first possibility of directly 
observing the general relativistic effects approaching the black hole horizon. 


CONTRIBUTION OF THE SHORT GRBS TO THE BLACK HOLE THEORY 


On the gravitational binding energy of white dwarf and neutron stars 


The aim of this section is to point out how the knowledge obtained from the black hole model is of relevance also 
for the basic theory of black holes and further how very high precision verification of general relativistic effects in the 
very strong field near the formation of the horizon should be expected in the near future. 

We shall first see how Eq.(63) for Mi„, 

Mto=Mo-|| + r+, (78) 

leads to a deeper physical understanding of the role of the gravitational interaction in the maximum energy extraction 
process of a black hole. This formula can also be of assistance in clarifying some long lasting epistemological issue 
on the role of general relativity, quantum theory and thermodynamics. 

It is well known that if a spherically symmetric mass distribution without any electromagnetic structure undergoes 
free gravitational collapse, its total mass-energy M is conserved according to the Birkhoff theorem: the increase in the 
kinetic energy of implosion is balanced by the increase in the gravitational energy of the system. If one considers the 
possibility that part of the kinetic energy of implosion is extracted then the situation is very different: configurations 
of smaller mass-energy and greater density can be attained without violating Birkhoff theorem. 

We illustrate our considerations with two examples: one has found confirmation from astrophysical observations, 
the other promises to be of relevance for gamma ray bursts (GRBs) (see [51]). Concerning the first example, it is 
well known from the work of [111] that at the endpoint of thermonuclear evolution, the gravitational collapse of a 
spherically symmetric star can be stopped by the Fermi pressure of the degenerate electron gas (white dwarf). A 
configuration of equilibrium can be found all the way up to the critical number of particles 


Acri, = 0.775- 


(79) 



where the factor 0.775 comes from the coefficient of the solution of the Lane-Emden equation with polytropic 

index n = 3, and nipi = is the Planck mass, niQ is the nucleon mass and /r the average number of electrons per 
nucleon. As the kinetic energy of implosion is carried away by radiation the star settles down to a configuration of 
mass 

M = Acritmo-t/, (80) 

where the gravitational binding energy U can be as high as 5.72 x lO^^AcritWo- 

Similarly Gamov (see [112]) has shown that a gravitational collapse process to still higher densities can be stopped 
by the Fermi pressure of the neutrons (neutron star) and Oppenheimer [113] has shown that, if the effects of strong 
interactions are neglected, a configuration of equilibrium exists also in this case all the way up to a critical number of 
particles 

Acrit = 0.398^, (81) 

“o 

where the factor 0.398 comes now from the integration of the Tolman-Oppenheimer-Volkoff equation (see e.g. [114]). 
If the kinetic energy of implosion is again carried away by radiation of photons or neutrinos and antineutrinos the final 
configuration is characterized by the formula (80) with U < 2.48 x These considerations and the existence 

of such large values of the gravitational binding energy have been at the heart of the explanation of astrophysical 
phenomena such as red-giant stars and supernovae: the corresponding measurements of the masses of neutron stars 
and white dwarfs have been carried out with unprecedented accuracy in binary systems [115]. 


On the minimum value of the reducible mass of a hlack hole formed in a spherically 

symmetric gravitational collapse 

From a theoretical physics point of view it is still an open question how far such a sequence can go: using causality 
nonviolating interactions, can one find a sequence of braking and energy extraction processes by which the density and 
the gravitational binding energy can increase indefinitely and the mass-energy of the collapsed object be reduced at 
will? This question can also be formulated in the mass-formula language of a black hole given in [11] (see also [51]): 
given a collapsing core of nucleons with a given rest mass-energy Mq, what is the minimum irreducible mass of the 
black hole which is formed? 

Following [50] and [51], consider a spherical shell of rest mass Mq collapsing in a flat space-time. In the neutral 
case the irreducible mass of the final black hole satisfies the equation (see [51]) 

Mi„=M = Mo-^ + r+, (82) 

where M is the total energy of the collapsing shell and T+ the kinetic energy at the horizon r+. Recall that the area S 
of the horizon is [11] 

S = Anr\ = 167rM^ (83) 

where r+ = 2Min- is the horizon radius. The minimum irreducible mass is obtained when the kinetic energy at 

the horizon T+ is 0, that is when the entire kinetic energy T+ has been extracted. We then obtain the simple result 

(84) 

We conclude that in the gravitational collapse of a spherical shell of rest mass Mq at rest at infinity (initial energy 
Mi = Mq), an energy up to 50% of Mqc^ can in principle be extracted, by braking processes of the kinetic energy. In 
this limiting case the shell crosses the horizon with T+ — 0. The limit ^ in the extractable kinetic energy can further 
increase if the collapsing shell is endowed with kinetic energy at infinity, since all that kinetic energy is in principle 
extractable. 

In order to illustrate the physical reasons for this result, using the formulas of [50], we have represented in Fig. 34 
the world lines of spherical shells of the same rest mass Mq, starting their gravitational collapse at rest at selected radii 
R*. These initial conditions can be implemented by performing suitable braking of the collapsing shell and concurrent 
kinetic energy extraction processes at progressively smaller radii (see also Fig. 35). The reason for the existence of 
the minimum (84) in the black hole mass is the “self closure” occurring by the formation of a horizon in the initial 
configuration (thick line in Fig. 34). 
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Figure 34. Collapse curves for neutral shells with rest mass Mq starting at rest at selected radii R* computed by using the exact 
solutions given in [50]. A different value of (and therefore of r+) corresponds to each curve. The time parameter is the 
Schwarzschild time coordinate t and the asymptotic behavior at the respective horizons is evident. The limiting configuration 
M[n- = ^ (solid line) corresponds to the case in which the shell is trapped, at the very beginning of its motion, by the formation of 
the horizon. 


Is the limit Min- ^ ^ actually attainable without violating causality? Let us consider a collapsing shell with charge 
Q. IfM > 2 a black hole is formed. As pointed out in [51] the irreducible mass of the final black hole does not depend 
on the charge Q. Therefore Eqs.(82) and (84) still hold in the charged case with r+ = M + a/M^ — Q^. In Fig. 35 
we consider the special case in which the shell is initially at rest at infinity, i.e. has initial energy M\ = Mq, for three 
different values of the charge Q. We plot the initial energy M,, the energy of the system when all the kinetic energy of 

implosion has been extracted as well as the sum of the rest mass energy and the gravitational binding energy —of 
the system (here R is the radius of the shell). In the extreme case Q = Mq, the shell is in equilibrium at all radii (see 
[50]) and the kinetic energy is identically zero. In all three cases, the sum of the extractable kinetic energy T and the 
electromagnetic energy ^ reaches 50% of the rest mass energy at the horizon, according to Eq.(84). 

What is the role of the electromagnetic field here? If we consider the case of a charged shell with Q ~ Mq, 
the electromagnetic repulsion implements the braking process and the extractable energy is entirely stored in the 
electromagnetic field surrounding the black hole (see [51]). In [51] we have outlined two different processes of 
electromagnetic energy extraction. We emphasize here that the extraction of 50% of the mass-energy of a black hole 
is not specifically linked to the electromagnetic field but depends on three factors: a) the increase of the gravitational 
energy during the collapse, b) the formation of a horizon, c) the reduction of the kinetic energy of implosion. Such 
conditions are naturally met during the formation of an extreme black hole but are more general and can indeed occur 
in a variety of different situations, e.g. during the formation of a Schwarzschild black hole by a suitable extraction of 
the kinetic energy of implosion (see Fig. 34 and Fig. 35). 


On the Bekenstein-Hawking consideration of incompatibility between general relativity and 

thermodynamics 

Now consider a test particle of mass m in the gravitational field of an already formed Schwarzschild black hole 
of mass M and go through such a sequence of braking and energy extraction processes. Kaplan [116] found for the 
energy E of the particle as a function of the radius r 


1-^ 


(85) 
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Figure 35. Energetics of a shell such that M; = Mq, for selected values of the charge. In the first diagram 2 = 0; the dashed line 
represents the total energy for a gravitational collapse without any braking process as a function of the radius R of the shell; the 
solid, stepwise line represents a collapse with suitable braking of the kinetic energy of implosion at selected radii; the dotted line 
represents the rest mass energy plus the gravitational binding energy. In the second and third diagram Q/Mq = 0.7, Q/Mq = I 
respectively; the dashed and the dotted lines have the same meaning as above; the solid lines represent the total energy minus the 
kinetic energy. The region between the solid line and the dotted line corresponds to the stored electromagnetic energy. The region 
between the dashed line and the solid line corresponds to the kinetic energy of collapse. In all the cases the sum of the kinetic energy 
and the electromagnetic energy at the horizon is 50% of Mq. Both the electromagnetic and the kinetic energy are extractable. It is 
most remarkable that the same underlying process occurs in the three cases: the role of the electromagnetic interaction is twofold: 
a) to reduce the kinetic energy of implosion by the Coulomb repulsion of the shell; b) to store such an energy in the region around 
the black hole. The stored electromagnetic energy is extractable as shown in [51]. 


It would appear from this formula that the entire energy of a particle could be extracted in the limit r 2M. Such 
100% efficiency of energy extraction has often been quoted as evidence for incompatibility between General Relativity 
and the second principle of Thermodynamics (see [117] and references therein). J. Bekenstein and S. Hawking have 
gone as far as to consider General Relativity not to be a complete theory and to conclude that in order to avoid 
inconsistencies with thermodynamics, the theory should be implemented through a quantum description [117, 118]. 
Einstein himself often expressed the opposite point of view (see e.g. [119]). 

The analytic treatment presented in [50] can clarify this fundamental issue. It allows to express the energy increase 
£ of a black hole of mass Mi through the accretion of a shell of mass Mq starting its motion at rest at a radius R in the 
following formula which generalizes Eq.(85): 


E=M-Mi 


( 86 ) 


where M — Mi + E is clearly the mass-energy of the final black hole. This formula differs from the Kaplan formula 
(85) in three respects: a) it takes into account the increase of the horizon area due to the accretion of the shell; b) it 
shows the role of the gravitational self energy of the imploding shell; c) it expresses the combined effects of a) and b) 
in an exact closed formula. 














The minimum value of E is attained for the minimum value of the radius R = 2M: the horizon of the final black 
hole. This corresponds to the maximum efficiency of the energy extraction. We have 
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(87) 


or solving the quadratic equation and choosing the positive solution for physical reasons 

Eram = \[\jM\+Ml-M^. 

The corresponding efficiency of energy extraction is 


^max 



( 88 ) 


(89) 


which is strictly smaller than 100% for any given Mq ^ 0. It is interesting that this analytic formula, in the limit 
M\ ^ Mo, properly reproduces the result of equation (84), corresponding to an efficiency of 50%. In the opposite limit 
M\ » Mo we have 

(90) 

Only for Mo —> 0, Eq.(89) corresponds to an efficiency of 100% and correctly represents the limiting reversible 
transformations introduced in [11]. It seems that the difficulties of reconciling General Relativity and Thermodynamics 
are ascribable not to an incompleteness of General Relativity but to the use of the Kaplan formula in a regime in which 
it is not valid. The generalization of the above results to stationary black holes is being considered. 


ON A SEPARATRIX IN THE GRAVITATIONAL COLLAPSE TO AN OVERCRITICAL 

ELECTROMAGNETIC BLACK HOLE 

We are now ready to analyze the dynamical properties of an electron-positron-photon plasma created by the vacuum 
polarization process occurring around a charged gravitationally collapsing core of an initially neutral star are examined 
within the framework of General Relativity and Quantum Field Theory. The Reissner-Nordstrom geometry is assumed 
to apply between the collapsing core and the oppositely charged remnant of the star. The appearance of a separatrix 
at radius R, well outside the asymptotic approach to the horizon, is evidenced. The neutral electron-positron-photon 
plasma created at radii r > R self-propels outwards to infinity, following the classical PEM-pulse analysis [55, 59]. 
The plasma created at r < R remains trapped and follows the gravitational collapse of the core only contributing to the 
reduction of the electromagnetic energy of the black hole and to the increase of its irreducible mass. This phenomenon 
has consequences for the observational properties of gamma-ray bursts and is especially relevant for the theoretical 
prediction of the temporal and spectral structure of the short bursts. 

The formulation of the physics of the dyadosphere of an electromagnetic black hole (black hole) has been until now 
approached by assuming the vacuum polarization process a la Sauter-Heisenberg-Euler-Schwinger [120, 25, 26] in 
the field of an already formed Kerr-Newmann [23] or Reissner-Nordstrom black hole [48, 51]. This acausal approach 
is certainly valid in order to describe the overall energetics and the time development of the gamma-ray bursts (GRBs) 
reaching a remarkable agreement between the observations and the theoretical prediction, in particular with respect to: 
a) the existence of a proper gamma-ray burst (P-GRB) [13], b) the afterglow detailed luminosity function and spectral 
properties [121, 28, 76] and c) the relative intensity of the P-GRB to the afterglow [14, 121, 28]. 

This acausal approach has to be improved by taking into account the causal dynamical process of the formation of 
the dyadosphere as soon as the detailed description on timescales of 10^^ — lO^^s of the P-GRB are considered. Such 
a description leads to theoretical predictions on the time variability of the P-GRB spectra which may become soon 
testable by a new class of specially conceived space missions. 

We report progress in this theoretically challenging process which is marked by distinctive and precise quantum 
and general relativistic effects. These new results have been made possible by the recent progress in Refs. [50], 
[51] and especially [104]. There it was demonstrated the intrinsic stability of the gravitational amplification of the 
electromagnetic field at the surface of a charged star core collapsing to a black hole. The plasma generated by 
the vacuum polarization process around the core is entangled in the electromagnetic field [49]. The pairs do 



thermalize in an electron-positron-photon plasma on a time scale 10^ — 10“^ times larger than h/nieC [104], where 
c is the speed of light and nig the electron mass. As soon as the thermalization has occurred, a dynamical phase of 
this electrically neutral plasma starts following the considerations already discussed in [55, 59]. While the temporal 
evolution of the 7 plasma takes place, the gravitationally collapsing core moves inwards, giving rise to a further 
amplified supercritical field, which in turn generates a larger amount of pairs leading to a yet higher temperature 
in the newly formed yplasma. We report, in the following, progress in the understanding of this crucial dynamical 
process: the main difference from the previous treatments is the fact that we do not consider an already formed black 
hole but we follow the dynamical phase of the formation of dyadosphere and of the asymptotic approach to the horizon 
by examining the time varying process at the surface of the gravitationally collapsing core. 

The space-time external to the surface of the spherically symmetric collapsing core is described by the Reissner- 
Nordstrom geometry [122] with line element 

ds^ = -a^dt^ + a^^dr^ + r^dD} , (91) 

with = dO^ + sin^ 9d<j)^, = a^{r) = I — 2M/r + jr^, where M and Q are the total energy and charge of 

the core as measured at infinity. On the core surface, which at the time to has radial coordinate ro, the electromagnetic 
field strength is (a = rf (ro) = G/r^. The equation of core’s collapse is (see [50]): 




(92) 


where H {yq) = ^ — and Mq is the core rest mass. Analytic expressions for the solution of Eq.(92) were given 

in [50]. We here recall that the dyadosphere radius is defined by S' (r^s) = Sc = jeh [48] as r^s = a/ eQh/m^S’, 
where e is the electron charge. In the following we assume that the dyadosphere starts to be formed at the instant 
^ds ~ Iq (Os) ~ 0- 

Having formulated the core collapse in General Relativity in Eq.(92), in order to describe the quantum phenomena, 
we consider, at each value of ro and to, a slab of constant coordinate thickness Ar small in comparison with r^s and 
larger than h/nieC^. All the results will be shown to be independent on the choice of the value of Ar. In each slab 
the process of vacuum polarization leading to pair creation is considered. As shown in [104, 49] the pairs 
created oscillate [123, 124, 125, 126] with ultrarelativistic velocities and partially annihilate into photons; the electric 
field oscillates around zero and the amplitude of such oscillations decreases with a characteristic time of the order of 
10^ — 10^ h/nieC^. The electric field is effectively screened to the critical value Sc and the pairs thermalize to an e+e ^7 
plasma. While the average of the electric field S over one oscillation is 0, the average of S'^ is of the order of S^, 
therefore the energy density in the pairs and photons, as a function of ro, is given by [51] 




-1 


(93) 


Eor the number densities of e+e pairs and photons at thermal equilibrium we have rig+g- ~ ny, correspondingly the 
equilibrium temperature Tq, which is clearly a function of yq and is different for each slab, is such that 


e (To) = fiy (7b) + EgY (7b) + Eg- {To) — eo, 


(94) 


with E and n given by Eermi (Bose) integrals (with zero chemical potential): 

(95) 

(96) 

where k is the Boltzmann constant. Erom the conditions set by Eqs.(94), (95), (96), we can now turn to the dynamical 
evolution of the plasma in each slab. We use the covariant conservation of energy momentum and the rate 

equation for the number of pairs in the Reissner-Nordstrom geometry external to the star core: 


(7b) = ^ 


e^{To) = ^{kTor 


(7b) = ^ [ ^^^,EdE, ny{To) = ^ {kTof 




= 0, Va {tlg+g-u'^) = av [nl+^- (r) - , 


(97) 



where = (e + p)u‘‘u^ + is the energy-momentum tensor of the plasma with proper energy density e and 
proper pressure p, is the fluid 4—velocity, n^+e- is the number of pairs, Wg+g- (T) is the equilibrium number of pairs 
and Wv is the mean of the product of the annihilation cross-section and the thermal velocity of pairs. We follow 
closely the treatment which we developed for the consideration of a plasma generated in the dyadosphere of an already 
formed black hole [55, 59]. It was shown in [55, 59] that the plasma expands as a pair-electromagnetic pulse (PEM 
pulse) of constant thickness in the laboratory frame. Since the expansion, hydrodynamical timescale is much larger 
than the pair creation {h/nieC^) and the thermalization ( 10 ^ — 10 ^/i/mgC^) time-scales, in each slab the plasma remains 
at thermal equilibrium in the initial phase of the expansion and the right hand side of the rate Eq.(97) is effectively 0, 
see Eig. 24 (second panel) of [28] for details. 

If we denote by the static Killing vector field normalized at unity at spacial infinity and by {E,}, the family 
of space-like hypersurfaces orthogonal to (f being the Killing time) in the Reissner-Nordstrom geometry, from 
Eqs.(97), the following integral conservation laws can be derived (see for instance [127, 128]) 

/ ^aT‘^^dI.h=E, [ n,+,-u'’dI.b=N,+,-, (98) 

Jz, JZ, 

where sin 9drd9d^ is the vector surface element, E the total energy and A(,+g- the total number of 

pairs which remain constant in each slab. We then have 

[{e + p)Y^-p]r^ = €, Hg+g-ya^V^ ^OTg+g-, (99) 

where € and 91g+g- are constants and 


7= a 


1 - a 


-4 



-1/2 


( 100 ) 


is the Lorentz 7 factor of the slab as measured by static observers. We can rewrite Eqs.(98) for each slab as 

^2 


(f) =«Vro, 

fWnf/gf pfrV 
[roj V®o/V"e+g-/ V“o/ eo[roJ 


( 101 ) 

( 102 ) 

(103) 


where pedex 0 refers to quantities evaluated at selected initial times to > 0 , having assumed r(to) = rg, drjdt = 0 , 
T (to) = To. 

Eq.(lOl) is only meaningful when /gg (r) > 0. Erom the structural analysis of such equation it is clearly identifiable 
a critical radius R such that: 


• for any slab initially located at ro > R we have /gg (r) > 0 for any value of r > tq and /gg (r) < 0 for r < tq; 
therefore a slab initially located at a radial coordinate ro > R moves outwards, 

• for any slab initially located at ro < R we have /gg (r) > 0 for any value of r+ < r < ro and /gg (r) < 0 for r > ro; 
therefore a slab initially located at a radial coordinate ro < R moves inwards and is trapped by the gravitational 
field of the collapsing core. 

We define the surface r = R, the dyadosphere trapping surface (DTS). The radius R of DTS is generally evaluated 
= 0. R is so close to the horizon value r+ that the initial temperature Tq satisfies kTo /$> nieC^ 

r=R 

and we can obtain for R an analytical expression. Namely the ultrarelativistic approximation of all Eermi integrals, 
Eqs.(95) and (96), is justified and we have Wg+g- {T) and therefore /gg ~ 1 — {T/Tof’ (oo/ct)^ {r < R). 

The defining equation of R, together with (103), then gives 


by the condition 


dr 


R^2M 


1 -f (1 - 2>Q^/ANf-) 


1/2 


>r+. 


(104) 


In the case of a black hole with M = 20Mq, Q = O.IM, we compute: 

• the fraction of energy trapped in DTS: 

E = f aeoc/E~0.53 f aeoc/E; 

Jrj^<r<R Jr+<r<r^s 


(105) 



• the world-lines of slabs of plasma for selected ro in the interval (see Fig. 36); 

• the world-lines of slabs of plasma for selected ro in the interval {r+,R) (see Fig. 37). 

At time T=to (R) when the DTS is formed, the plasma extends over a region of space which is almost one order of 
magnitude larger than the dyadosphere and which we define as the effective dy ado sphere. The values of the Lorentz 7 
factor, the temperature and number density in the effective dyadosphere are given in Fig. 38. 

In conclusion we see how the causal description of the dyadosphere formation can carry important messages on the 
time variability and spectral distribution of the P-GRB due to quantum effects as well as precise signature of General 
Relativity. 


DESCRIPTION OF THE ELECTRON-POSITRON PLASMA OSCILLATIONS BY 
GENERALIZED VLASOV-BOLTZMANN-MAXWELL EQUATION IN THE PEM 

PULSE PHASE 

We describe the creation and evolution of electron-positron pairs in a strong electric field as well as the pairs 
annihilation into photons. The formalism is based on generalized Vlasov equations, which are numerically integrated. 
We recover previous results about the oscillations of the charges, discuss the electric field screening and the relaxation 
of the system to a thermal equilibrium configuration. The timescale of the thermalization is estimated to be ~ 
10 ^ — 10 “^ hlnieC^. 

On the observability of electron-positron pairs created in vacuum polarization in Earth 

bound experiment and in astrophysics 

Three different earth-bound experiments and one astrophysical observation have been proposed for identifying the 
polarization of the electronic vacuum due to a supercritical electric field (^ > jeh, where nie and e are the 

electron mass and charge) postulated by Sauter-Heisenberg-Euler-Schwinger [120]: 

1. In central collisions of heavy ions near the Coulomb barrier, as first proposed in [129, 130] (see also [131, 
132, 133]). Despite some apparently encouraging results [134], such efforts have failed so far due to the small 
contact time of the colliding ions [135, 136, 137, 138, 139]. Typically the electromagnetic energy involved in the 
collisions of heavy ions with impact parameter h ^ lO^^^cm is £1 ~ lO^^erg and the lifetime of the diatomic 
system is t\ ~ lO^^^s. 

2. In collisions of an electron beam with optical laser pulses: a signal of positrons above background has been 
observed in collisions of a 46.6 GeV electron beam with terawatt pulses of optical laser in an experiment at the 
Final Focus Test Beam at SLAC [140]; it is not clear if this experimental result is an evidence for the vacuum 
polarization phenomenon. The energy of the laser pulses was £2 ^ lO^erg, concentrated in a space-time region 
of spacial linear extension (focal length) I 2 ~ lO^^cm and temporal extension (pulse duration) t 2 ~ 10^'^s [140]. 

3. At the focus of an X-ray free electron laser (XFEL) (see [141, 142, 143] and references therein). Proposals for 
this experiment exist at the TESLA collider at DESY and at the LCLS facility at SLAC [141]. Typically the 
electromagnetic energy at the focus of an XEEL can be £3 ~ 10®erg, concentrated in a space-time region of 
spacial linear extension (spot radius) Z 3 ~ lO^^cm and temporal extension (coherent spike length) t^ ~ 10 ^'^s 
[141]. 

and from astrophysics 

1. around an electromagnetic black hole (black hole) [23, 48, 144], giving rise to the observed phenomenon of 
gamma-ray bursts (GRB) [13, 14, 15, 76]. The electromagnetic energy of an black hole of mass M ^ IQMq 
and charge 2 ~ 0 . IM/\fG is £4 ^ 10 ^‘^erg and it is deposited in a space-time region of spacial linear extension 
I 4 ~ lO^cm [48, 51] and temporal extension (collapse time) f 4 ~ lO^^s [104]. 




Figure 36. World line of the collapsing charged core (dashed line) as derived from Eq.(92) for a hlack hole with M = 20Mq, 
Q = 0. IM; world lines of slabs of plasma for selected radii rg in the interval (R,rg^). At time t the expanding plasma extends over 
a region which is almost one order of magnitude larger than the dyadosphere. The small rectangle in the right bottom is enlarged in 
Fig. 37. 


On the role of transparency condition in the electron-positron plasma 

In addition to their marked quantitative difference in testing the same basic physical phenomenon, there is a very 
important conceptual difference among these processes: the first three occur in a transparency condition in which 
the created electron-positron pairs and, possibly, photons freely propagate to infinity, while the one in the black hole 










Figure 37. Enlargement of the small rectangle in the right bottom of Fig. 36. World-lines of slabs of plasma for selected radii rg 
in the interval (r+,.R). 


occurs in an opacity condition [59]. Under the opacity condition a thermalization effect occurs and a final equipartition 
between the and 7 is reached. Far from being just an academic issue, this process and its characteristic timescale is 
of the greatest importance in physics and astrophysics. It has been shown by a numerical simulation done in Livermore 
and an analytic work done in Rome [59], that, as soon as the thermalization of and / created around a black hole 
has been reached, the plasma self propels outwards and this process is at the very heart of the gamma-ray burst (GRB) 
phenomenon. A critical step was missing up to now: how to bridge the gap between the creation of pairs in the 
supercritical field of the black hole and the thermalization of the system to a plasma configuration. We report some 
progress on this topic with special attention to the timescale needed for the thermalization of the newly created 
pairs in the background field. The comparison of the thermalization timescale to the one of gravitational collapse, 
which occurs on general relativistic timescale, is at the very ground of the comprehension of GRBs [104]. 

The evolution of a system of particle-antiparticle pairs created by the Schwinger process has been often described by 
a transport Vlasov equation (see, for example, [145, 146]). More recently it has been showed that such an equation can 
be derived from quantum field theory [147, 148, 149]. In the homogeneous case, the equations have been numerically 
integrated taking into account the back reaction on the external electric field [123, 124, 125, 126]. In many papers 
(see [150] and references therein) a phenomenological term describing equilibrating collisions is introduced in the 
transport equation which is parameterized by an effective relaxation time T. In [150] one further step is taken by 
allowing time variability of t; the ignorance on the collision term is then parameterized by a free dimensionless 
constant. The introduction of a relaxation time corresponds to the assumption that the system rapidly evolves towards 
thermal equilibrium. In this paper we focus on the evolution of a system of pairs, explicitly taking into account 
the scattering processes ^ 77 . Since we are mainly interested in a system in which the electric field varies on 
macroscopic length scale (I ~ lO^cm, above), we can limit ourselves to a homogeneous electric field. Also, we will use 
transport equations for electrons, positrons and photons, with collision terms, coupled to Maxwell equations. There is 
no free parameter here: the collision terms can be exactly computed, since the QED cross sections are known. Starting 
from a regime which is far from thermal equilibrium, we find that collisions do not prevent plasma oscillations in the 
initial phase of the evolution and analyse the issue of the timescale of the approach to a 7 plasma equilibrium 
configuration, which is the most relevant quantity in the process of gravitational collapse [104]. 


The Vlasov-Boltzmann-Maxwell equations and their solutions 

The motion of positrons (electrons) is the resultant of three contributions: the pair creation, the electric acceleration 
and the annihilation damping. The homogeneous system consisting of electric field, electrons, positrons and photons 














Figure 38. Physical parameters in the effective dyadosphere: Lorentz y factor, proper temperature and proper e+e number 
density as functions at time t for a black hole with M = 20Mq and Q = 0. IM. 


can be described by the equations 


dtfe + eEd^fe = y (E,p) - (f,p), 

(106) 


(107) 

5,E = -jp(E)-jc(f), 

(108) 


where = /e(f,p) is the distribution function in the phase-space of positrons (electrons), fy = /y(t,k) is the 
distribution function in the phase-space of photons, E is the electric field, Ep = (p • p + ni^') ^ is the energy of an 
electron of 3-momentum p {me is the mass of the electron) and et = (k • k) is the energy of a photon of 3-momentum 
k. fe and fy are normalized so that / fe (t,p) = (f), / /y (f,k) = ny{t) , where rig and riy are number 

(ZTT) 




















densities of positrons (electrons) and photons, respectively. The term 

= (2^)' = - kEllog [l -exp(-^^(^)] 5(p||) (109) 

is the Schwinger source for pair creation (see [123, 124]): p|| and px are the components of the 3-momentum p 
parallel and orthogonal to E. We assume that the pairs are produced at rest in the direction parallel to the electric field 
[123, 124], We also have, in Eqs. (106), (107) and (108), 




X [/«(p)/<"(pi)-/r(ki)/r(k2)], 

(110) 



X [/e(Pl)/<"(P2)-/y(k)/y(ki)] , 

(111) 


which describe probability rates for pair creation by photons and pair annihilation into photons, ^ = 
•^c+(Pi)e-(P 2 )^y{k) 7 (ki) being the matrix element for the process e^(pi)e^(p 2 ) ^ 7 (k) 7 (ki). Note that the 
collisional terms (110) and (111) are either inapplicable or negligible in the case of the above three earth-bound 
experiments where the created pairs do not originate a dense plasma. They have been correctly neglected in previous 
works (see e. g. [143]). Collisional terms have also been considered in the different physical context of vacuum 
polarization by strong chromoelectric fields. Unlike the present QED case, where expressions for the cross sections 
are known exactly, in the QCD case the cross sections are yet unknown and such collisional terms are of a phe¬ 
nomenological type and useful uniquely near the equilibrium regime [150]. Einally jp (E) = 

and jc(f) = 2eneJ ^ ^3 j-fe (p) are polarization and conduction current respectively (see [146]). In Eqs. (110) 

and (111) we neglect, as a first approximation, Pauli blocking and Bose enhancement (see e.g. [124]). By suitably 
integrating (106) and (107) over the phase spaces of positrons (electrons) and photons, we find the following exact 
equations for mean values: 

j^rie = 5(E) - nl (cTiv')^ -f <|(72v")y. 

(Cp)^ = erieE ■ (v)^ -f ^E • jp - (epdi v")^ + (ek(T2v")^, 

inr{£k)y = 2nl {e^Oiv')^-2n]{ei,02v") 


(p)g = eneE - n] (pcJi v')^ , 





|■E = -2eng(v)^-jp(E), 




( 112 ) 

where, for any function of the momenta 





(^<p' . 

^(pi,- 

•■.P«)-/«(Pi)-- 

fe (Pn) , 

(113) 

{C(k„...,k,)), = x7l5^i-W 

G(ki,. 

•■>k/)-/y(ki)-. 

•••/y(k/). 

(114) 


Eurthermore v' is the relative velocity between electrons and positrons, v" is the relative velocity between photons, 
(ji = ffi (gp is the total cross section for the process ^ 77 and 02 — 02 (gf cross section for 

the process 77 ^ (here gCoM energy of a particle in the reference frame of the center of mass). 

In order to evaluate the mean values in system (112) we need some further hypotheses on the distribution functions. 
Let us define p||, gp and such that {p\\)^ = p\\, (£p)g = gp = (p| + Px + We assume 

fe (f ,p) rie (t) 5 (p|| - P||) 5 (pi - pi). (115) 

Since in the scattering —> 77 the coincidence of the scattering direction with the incidence direction is statistically 

favored, we also assume 

/y(f,k) oc ny{t) 5 (ki -ki) [5 (k|| -k||) -f 5 (k|| +k||)] , 


(116) 



where A:|| and kx have analogous meaning as ;?|| and px and the terms 5 (A:|| — ^||) and 5 (A:|| +k\\) account for the 
probability of producing, respectively, forwardly scattered and backwardly scattered photons. Since the Schwinger 
source term (109) implies that the positrons (electrons) have initially fixed P\\, P\\ = 0, assumption (115) ((116)) means 
that the distribution of p|| (A:||) does not spread too much with time and, analogously, that the distribution of energies is 
sufficiently peaked to be describable by a 5—function. The dependence on the momentum of the distribution functions 
has been discussed in [124, 148]. Approximations (115), (116) reduce Eqs. (112) to a system of ordinary differential 
equations. In average, since the inertial reference frame we fix coincides with the center of mass frame for the processes 
^ 77 , ~ g for each species. Substituting (115) and (116) into (112) we find 

= 5(^)-2n^(Jip7' \ 7tg\\\+2nya2, 

^ny = Anluip;^ I - 4n^a2, 

■jjPe = Pg I I “f 2^7/7 '2ngPgCT\Pg Itt^H |-|-2nyPyC72, 

^Py = AtlgPeOlPg I TTgll I — 4nyPyO'2, 

^2ng7rg|| (Jip^ I TT^ii I, 

^^ = -2enep;^\7t,\\\-jpi^), (117) 

where Pe = n^ep, Py — nyfik, T^e\\ — ”eP|| the energy density of positrons (electrons), the energy density of photons 
and the density of “parallel momentum” of positrons (electrons). S' is the electric field strength and jp the unique 
component of jp parallel to E. ai and 02 are evaluated at = g for each species. Note that Eqs.(117) are 
“classical” in the sense that the only quantum information is encoded in the terms describing pair creation and 
scattering probabilities. Eqs.(117) are consistent with energy density conservation: ^ (pg + Py + jS^) = 0. 

The initial conditions for Eqs.(117) are He = ny = Pe = Py = =Q, S = Sq. In Eig. 39 the results of the 

numerical integration for Sq = 9Sc is showed. The integration stops at f = 150 Tc (where Tc = h/rrigC^). Each 
variable is represented in units of nig and Ac = h/nigC. The numerical integration confirms [123, 124] that the system 
undergoes plasma oscillations: a) the electric field oscillates with decreasing amplitude rather than abruptly reaching 
the equilibrium value; b) electrons and positrons oscillates in the electric field direction, reaching ultrarelativistic 
velocities; c) the role of the ^ 77 scatterings is marginal in the early time of the evolution, the electrons are too 
extremely relativistic and consequently the density of photons builds up very slowly (see. details in Eig. 39). 

At late times the system is expected to relax to a plasma configuration of thermal equilibrium and assumptions 
(115) and (116) have to be generalized to take into account quantum spreading of the distribution functions. It is 
nevertheless interesting to look at the solutions of Eqs.(117) in this regime. In Eig. 40 we plot the numerical solution 
of Eqs.(l 17) but the integration extends here all the way up to f = 7000 Tc (the time scale of oscillations is not resolved 
in these plots). It is interesting that the leading term recovers the expected asymptotic behavior: a) the electric field is 
screened to about the critical value: S Sc for t ~ 10^ — lO'^Tc ^ Tc; b) the initial electromagnetic energy density is 
distributed over electron-positron pairs and photons, indicating energy equipartition; c) photons and electron-positron 
pairs number densities are asymptotically comparable, indicating number equipartition. At such late times a regime 
of thermalized electrons-positrons-photons plasma begins and the system is describable by hydrodynamic equations 
[104, 59]. 

We provided a very simple formalism apt to describe simultaneously the creation of electron-positron pairs by 
a strong electric field S > Sc and the pairs annihilation into photons. As discussed in literature, we find plasma 
oscillations. In particular the collisions do not prevent such a feature. This is because the momentum of electrons 
(positrons) is very high, therefore the cross section for the process ^ 77 is small and the annihilation into 
photons is negligible in the very first phase of the evolution. As a result, the system takes some time (f ~ 10^ — 10“^ Tc) to 
thermalize to a 7 plasma equilibrium configuration. We finally remark that, at least in the case of electromagnetic 

Schwinger mechanism, the picture could be quite different from the one previously depicted in literature, where the 
system is assumed to thermalize in a very short time (see [150] and references therein). 

It is conceivable that in the race to first identify the vacuum polarization process a la Sauter-Euler-Heisenberg- 
Schwinger, the astrophysical observations will reach a positive result before earth-bound experiments, much like in 
the case of the discovery of lines in the Sun chromosphere by J. N. Lockyer in 1869, later identified with the Helium 
spectral lines by W. Ramsay in 1895 [151]. 




Figure 39. Plasma oscillations. We set £0 = 9S’c, t < 150 tc and plot: a) electromagnetic field strength; b) electrons energy density; 
c) electrons number density; d) photons energy density; e) photons number density as functions of time. 


OBSERVATIONAL SIGNATURES OF AN ELECTROMAGNETIC OVERCRITICAL 
GRAVITATIONAL COLLAPSE AND PREDICTION OF SPECTRAL EVOLUTION OF 

SHORT GRBS 


We finally present theoretical predictions for the spectral, temporal and intensity signatures of the electromagnetic 
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Figure 40. Plasma oscillations. We set <Sq = 9<ac, t < 7000tc plot: a) electromagnetic field strength; b) electrons energy 
density; c) electrons number density; d) photons energy density; e) photons number density as functions of time - the oscillation 
period is not resolved in these plots. The model used should have a breakdown at a time much earlier than 7000Tc and therefore 
this plot contains no more than qualitative informations. 


radiation emitted during the process of the gravitational collapse of a stellar core to a black hole, during which 
electromagnetic field strengths rise over the critical value for pair creation. The last phases of this gravitational 
collapse are studied, using the result presented in the previous sections, leading to the formation of a black hole with 
a subcritical electromagnetic field, likely with zero charge, and an outgoing pulse of initially optically thick 
photon plasma. Such a pulse reaches transparency at Lorentz gamma factors of 10^-10^. We find a clear signature in 
the outgoing electromagnetic signal, drifting from a soft to a hard spectrum, on very precise time-scales and with a 







very specific intensity modulation. 

We outline finally the relevance of these theoretical results for the understanding of short gamma-ray bursts. 


The model 

The dynamics of the collapse of an electrically-charged stellar core, separating itself from an oppositely charged 
remnant in an initially neutral star, was first modeled by an exact solution of the Einstein-Maxwell equations corre¬ 
sponding to a shell of charged matter in Ref. [50]. The fundamental dynamical equations and their analytic solutions 
were obtained, revealing the amplification of the electromagnetic field strength during the process of collapse and 
the asymptotic approach to the final static configuration. The results, which properly account for general relativistic 
effects, are summarized in Fig. 1 and Fig. 2 of Ref. [50]. 

A first step toward the understanding of the process of extracting energy from a black hole was obtained in Ref. [51], 
where it was shown how the extractable electromagnetic energy is not stored behind the horizon but is actually 
distributed all around the black hole. Such a stored energy is in principle extractable, very efficiently, on time-scales 
~ by a vacuum polarization process a la Sauter-Heisenberg-Fuler-Schwinger [120, 25, 26]. Such a process 

occurs if the electromagnetic field becomes larger than the critical field strength for pair creation. In Ref. [51] 
we followed the approach of Damour and Ruffini [23] in order to evaluate the energy density and the temperature of 
the created -photon plasma. As a byproduct, a formula for the irreducible mass of a black hole was also derived 
solely in terms of the gravitational, kinetic and rest mass energies of the collapsing core. This surprising result allowed 
us in Ref. [52] to obtain a deeper understanding of the maximum limit for the extractable energy during the process 
of gravitational collapse, namely 50% of the initial energy of the star: the well known result of a 50% maximum 
efficiency for energy extraction in the case of a Reissner-Nordstrom black hole [11] then becomes a particular case of 
a process of much more general validity. 

The crucial issue of the survival of the electric charge of the collapsing core in the presence of a copious process 
of pair creation was addressed in Refs. [49, 53]. By using theoretical techniques borrowed from plasma physics 
and statistical mechanics [146, 123, 124, 125, 148, 149, 126] based on a generalized Vlasov equation, it was possible 
to show that while the core keeps collapsing, the created pairs are entangled in the overcritical electric field. The 
electric field itself, due to the back reaction of the created pairs, undergoes damped oscillations in sign finally 
settling down to the critical value ^c- The pairs fully thermalize to an e+e^-photon plasma on time-scales typically 
of the order of lO^-lO'^fi/nieC^. During this characteristic damping time, which we recall is much larger than the pair 
creation time-scale h/nieC^, the core moves inwards, collapsing with a speed 0.2-0.8c, further amplifying the electric 
field strength at its surface and enhancing the pair creation process. 

Turning now to the dynamical evolution of such an plasma we recall that, after some original attempt to 
consider a steady state emission [152, 153], the crucial progress was represented by the understanding that during the 
optically thick phase such a plasma expands as a thin shell. There exists a fundamental relation between the width 
of the expanding shell and the Lorentz gamma factor. The shell expands, but the Lorentz contraction is such that its 
width in laboratory frame appears to be constant. Such a result was found in [56] on the basis of a numerical approach, 
further analyzed in Bisnovatyi-Kogan and Murzina [57] on the basis of an analytic approach. Attention to the role of 
the rate equations governing the annihilation were given in [154], where approximations to the full equation were 
introduced. These results were improved in two important respects in 1999 and 2000 [55, 59]: the initial conditions 
were made more accurate by the considerations of the dyadosphere as well as the dynamics of the shell was improved 
by the self-consistent solution of the hydrodynamical equation and the rate equation for the plasma following 
both an analytic and numerical approach. 

We are now ready to report the result of using the approach in [55, 59] in this general framework describing the 
dynamical formation of the dyadosphere. 

The first attempt to analyze the expansion of the newly generated and thermalized -photon plasma was made 
in Ref. [49]. The initial dynamical phases of the expansion were analyzed, using the general relativistic equations of 
Ref. [50] for the gravitational collapse of the core. A separatrix was found in the motion of the plasma at a critical 
radius R\ the plasma created at radii larger than R expands to infinity, while the one created at radii smaller than R 
is trapped by the gravitational field of the collapsing core and implodes towards the black hole. The value of R was 
found in Ref. [49] to be R = 2GM/c^[\ + (l — where M and Q are the mass and the charge of the 

core, respectively. 

We now pursue further the evolution of such a system, describing the dynamical phase of the expansion of the pulse 



of the optically thick plasma all the way to the point where the transparency condition is reached. Some pioneering 
work in this respect were presented in Goodman in 1986 [155]. In this process the pulse reaches ultrarelativistic 
regimes with Lorentz factor 7 ~ 10^-10^. The spectra, the luminosities and the time-sequences of the electromagnetic 
signals captured by a far-away observer are analyzed here in detail for the first time. The relevance of these theoretical 
results for short-bursts is then discussed. 


The expansion of the e^e y plasma as a discrete set of elementary slabs 


We discretize the gravitational collapse of a spherically symmetric core of mass M and charge Q by considering a 
set of events along the world line of a point of fixed angular position on the collapsing core surface. Between each of 
these events we consider a spherical shell slab of plasma of constant coordinate thickness Ar so that: 

1. Ar is assumed to be a constant which is small with respect to the core radius; 

2. Ar is assumed to be large with respect to the mean free path of the particles so that the statistical description of 
the 7 plasma can be used; 

3. There is no overlap among the slabs and their union describes the entirety of the process. 

We check that the final results are independent of the special value of the chosen Ar. 

In order to describe the dynamics of the expanding plasma pulse the energy-momentum conservation law and the 
rate equation for the number of pairs in the Reissner-Nordstrom geometry external to the collapsing core have to be 
integrated: 


(«.+.-[«?+.- 


(118) 

(119) 


where T^'' = {e-\- p) u^u'' + is the energy-momentum tensor of the plasma with proper energy density e and 
proper pressure p, u“ is the fluid 4-velocity, n^+g- is the pair number density, n^+g- (T) is the equilibrium pair number 
density at the temperature T of the plasma and oT is the mean of the product of the annihilation cross-section 
and the thermal velocity of the pairs. We use Eqs.(118) and (119) to study the expansion of each slab, following 
closely the treatment developed in Refs [55, 59] where it was shown how a homogeneous slab of plasma expands as a 
pair-electromagnetic pulse (PEM pulse) of constant thickness in the laboratory frame. Two regimes can be identified 
in the expansion of the slabs: 


1. In the initial phase of expansion the plasma experiences the strong gravitational field of the core and a fully 
general relativistic description of its motion is needed. The plasma is sufficiently hot in this first phase that the 
pairs and the photons remain at thermal equilibrium in it. As shown in Ref. [49], under these circumstances, 
the right hand side of Eqs.(l 19) is effectively 0 and Eqs.(l 18) and (119) are equivalent to: 



( 120 ) 


1 

where r is the radial coordinate of a slab of plasma, a = (l — 2MG/c^r +Q^G/c^r^') ' is the gravitational 
redshift factor and the subscript “o" refers to quantities evaluated at the initial time. 

2. At asymptotically late times the temperature of the plasma drops below an equivalent energy of 0.5 MeV and 
the pairs and the photons can no longer be considered to be in equilibrium: the full rate equation for pair 
annihilation needs to be used. However, the plasma is so far from the central core that gravitational effects can be 
neglected. In this new regime, as shown in Ref. [55], Eqs.(l 18) and (119) reduce to: 
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Figure 41. Expansion of the plasma created around an overcritical collapsing stellar core with M = lOM© and Q = Q.I'/GM. 
Upper diagram: world lines of the plasma. Lower diagram: Lorentz y factor as a function of the radial coordinate r. 


where F = 1 + p/e, 'f' is the volume of a single slab as measured in the laboratory frame by an observer at rest 
with the black hole, Aig+g- = yUg+g- is the pair number density as measured in the laboratory frame by an observer 
at rest with the black hole, and Aig+g- (F) is the equilibrium laboratory pair number density. 


The reaching of transparency and the signature of the outgoing gamma ray signal 

Eqs.(120) and (121) must be separately integrated and the solutions matched at the transition between the two 
regimes. The integration stops when each slab of plasma reaches the optical transparency condition given by 

rAr 

/ (JT-ng+g-c/r ~ 1, (122) 

Jo 

where Gt is the Thomson cross-section and the integral extends over the radial thickness Ar of the slab. The evolution 
of each slab occurs without any collision or interaction with the other slabs; see the upper diagram in Fig. 41. The 
outer layers are colder than the inner ones and therefore reach transparency earlier; see the lower diagram in Fig. 41. 
In Fig. 41, Fqs.(120) and (121) have been integrated for a core with 


M=10Mq, Q = 0.iVgM; (123) 

the upper diagram represents the world lines of the plasma as functions of the radius, while the lower diagram shows 
the corresponding Lorentz y factors. The overall independence of the result of the dynamics on the number N of the 
slabs adopted in the discretization process or analogously on the value of Ar has also been checked. We have repeated 
the integration for A^ = 10, A^ = 100 reaching the same result to extremely good accuracy. The results in Fig. 41 
correspond to the case A^ = 10 . 

We now turn to the results in Fig. 42, where we plot both the theoretically predicted luminosity L and the spectral 
hardness of the signal reaching a far-away observer as functions of the arrival time ta- Since all three of these quantities 
depend in an essential way on the cosmological redshift factor z, see Refs. [156, 121], we have adopted a cosmological 
redshift z = 1 for this figure. 

As the plasma becomes transparent, gamma ray photons are emitted. The energy hco of the observed photon is 
h(0 = kyT / (1 -fz), where k is the Boltzmann constant, T is the temperature in the comoving frame of the pulse and 7 
is the Lorentz factor of the plasma at the transparency time. We also recall that if the initial zero of time is chosen as 
the time when the first photon is observed, then the arrival time ta of a photon at the detector in spherical coordinates 






Figure 42. Predicted observed luminosity and observed spectral hardness of the electromagnetic signal from the gravitational 
collapse of a collapsing core with M = IOMq, Q = Q.l\fGM at z = 1 as functions of the arrival time ta- 


centered on the black hole is given by [156, 121]: 

f« = (l+z)[f + ^-^cos0] (124) 

where {t ,r{t) labels the laboratory emission event along the world line of the emitting slab and rp is the initial 
position of the slab. The projection of the plot in Fig. 42 onto the ta-L plane gives the total luminosity as the sum of the 
partial luminosities of the single slabs. The sudden decrease of the intensity at the time t = 0.040466 s corresponds to 
the creation of the sepamtrix introduced in Ref. [53]. We find that the duration of the electromagnetic signal emitted 
by the relativistically expanding pulse is given in arrival time by 

Afa~5xl0^2s. (125) 

The projection of the plot in Fig. 42 onto the kTobs, ta plane describes the temporal evolution of the spectral hardness. 
We observe a precise soft-to-hard evolution of the spectrum of the gamma ray signal from ~ 10^ KeV monotonically 
increasing to ~ 1 MeV. We recall that kTobs = kjT / (1 +z). 

The above quantities are clearly functions of the cosmological redshift z, of the charge Q and the mass M of the 
collapsing core. We present in Fig. 3 the arrival time interval for M ranging from M ~ lOM© to lO^M©, keeping 
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Figure 43. Arrival time duration of the electromagnetic signal from the gravitational collapse of a stellar core with charge 
Q = Q.i\fGM as a function of the mass M of the core. 


2 = 0.1 ^/GM. The arrival time interval is very sensitive to the mass of the black hole: 

Ata~10^^-10^'s. (126) 

Similarly the spectral hardness of the signal is sensitive to the ratio Qj^fCM [54]. Moreover the duration, the spectral 
hardness and luminosity are all sensitive to the cosmological redshift z (see Ref. [54]). All the above quantities can 
also be sensitive to a possible baryonic contamination of the plasma due to the remnant of the progenitor star which 
has undergone the process of gravitational collapse. 


CONCLUSIONS 

On the GRB-Supernova connection 

We first stress some general considerations originating from comparing and contrasting the three GRB sources we 
have discussed: 

1. The value of the B parameter for all three sources occurs, as theoretically expected, in the allowed range (see Fig. 

22 ) 

10^^<B<10^^ (127) 

We have in fact: 

GRB 991216 B = 3.0x E^ya = 4.8 x 10^^ erg 

GRB 980425 B = 7.0 x E^ya = 1.1 x 10"^*^ erg 

GRB 030329 B = 4.8 x lO^^ E^ya = 2.1 x 10^2 g^g 

2. The enormous difference in the GRB energy of the sources simply relates to the electromagnetic energy of the 

black hole given in Eq.(6) which turns out to be smaller than the critical value given by Eq.(7). The fact that the 
theory is valid over 5 orders of magnitude is indeed very satisfactory. 

3. Also revealing is the fact that in both sources GRB 980425 and GRB 030329 the associated supernova energies 
are similar. We have, in fact, for both SN 1998bw and SN 2003dh an energy ~ 10^® erg. Details in Eraschetti 
et al. [17] and Bernardini et al. [97]. The further comparison between the SN luminosity and the GRB intensity 
is crucial. In the case of GRB 980425 the GRB and the SN energies are comparable, and no dominance of one 




source over the other can be ascertained. In the case of GRB 030329 the energy of the GRB source is 10^ larger 
than the SN: in no way the GRB can originate from the SN event. 

The above stringent energetics considerations and the fact that GRBs occur also without an observed supernova give a 
strong evidence that GRBs cannot originate from supernovae. 


URCA-1 and URCA-2 

We turn now to the most exciting search for the nature of URCA-1 and URCA-2. We have already mentioned above 
that a variety of possibilities naturally appear. The first possibility is that the URCA sources are related to the black 
hole originating the GRB phenomenon. In order to probe such an hypothesis, it would be very important to find even 
a single case in which an URCA source occurs in association with a GRB and in absence of an associated supernova. 
Such a result, theoretically unexpected, would open an entire new problematic in relativistic astrophysics and in the 
physics of black holes. 

If indeed, as we expect, the clear association between URCA sources and the supernovae occurring together with 
the GRBs, then it is clear that the analysis of the other two possibilities will be favored. Namely, an emission from 
processes occurring in the early phases of the expansion of the supernova remnant or the very exciting possibility 
that for the first time we are observing a newly born neutron star out of the supernova phenomenon. Of course, this 
last hypothesis is the most important one, since it would offer new fundamental information about the outcome of the 
gravitational collapse, about the equations of state at supranuclear densities and about a variety of fundamental issues 
of relativistic astrophysics of neutron stars. We shall focus in the following only on this last topic. 

We have already recalled how the need for a rapid cooling process due to neutrino anti-neutrino emission in the 
process of gravitational collapse leading to the formation of a neutron star was considered for the first time by George 
Gamow and Mario Schoenberg in 1941 [157]. It was Gamow who gave this process the name “Urea process”, see 
and . Since then, a systematic analysis of the theory of neutron star cooling was advanced by Tsuruta [158, 159], 
Tsuruta and Cameron [160], Tsuruta et al. [161] and by Canuto [162]. The coming of age of X-ray observatories such 
as Einstein (1978-1981), EXOSAT (1983-1986), ROSAT (1990-1998), and the contemporary missions of Chandra 
and XMM-Newton since 1999 dramatically presented an observational situation establishing very embarrassing and 
stringent upper limits to the surface temperature of neutron stars in well known historical supernova remnants (see e.g. 
Romani [163]). It was so that, for some remnants, notably SN 1006 and the Tycho supernova, the upper limits to the 
surface temperatures were significantly lower than the temperatures given by standard cooling times (see e.g. Romani 
[163]). Much of the theoretical works has been mainly directed, therefore, to find theoretical arguments in order to 
explain such low surface temperature 7^ ~ 0.5-1.0 x 10® K — embarrassingly low, when compared to the initial hot 
(~ 10^* K) birth of a neutron star in a supernova explosion (see e.g. Romani [163]). Some very important steps in this 
direction of research have been represented by the works of Van Riper [164, 165], Lattimer and his group [166, 167] 
and by the most extensive work of Yakovlev and his group [168]. The youngest neutron star to be searched for using 
its thermal emission in this context has been the pulsar PSR J0205H-6449 in 3C 58 (see e.g. Yakovlev and Pethick 
[168]), which is 820 years old! Recently, evidence for the detection of thermal emission from the crab nebula pulsar 
was reported by Trumper [169] which is, again, 951 years old. 

In the case of URCA-1 and URCA-2, we are exploring a totally different regime: the X-ray emission possibly from 
a recently born neutron star in the first days - months of its existence, where no observations have yet been performed 
and no embarrassing constraints upper limits on the surface temperature exist. The reason of approaching first the issue 
of the thermal emission from the neutron star surface is extremely important, since in principle it can give information 
on the equations of state in the core at supranuclear densities and on the detailed mechanism of the formation of the 
neutron star itself and the related neutrino emission. It is of course possible that the neutron star is initially fast rotating 
and its early emission is dominated by the magnetospheric emission or by accretion processes from the remnant which 
would overshadow the thermal emission. In that case a periodic signal related to the neutron star rotational period 
should in principle be observable in a close enough GRB source provided the suitable instrumentation from the Earth. 

The literature on young born neutron star is relatively scarse today. There are some very interesting contributions 
which state: “The time for a neutron star’s center to cool by the direct URCA process to a temperature T has been 
estimated to be f = 20 \T / (lO^X)] ^ s. The direct URCA process and all the exotic cooling mechanisms only occur 
at supranuclear densities. Matter at subnuclear densities in neutron star crust cools primarily by diffusion of heat to the 
interior. Thus the surface temperature remains high, in the vicinity of 10® K or more, until the crust’s heat reservoir is 
consumed. After this diffusion time, which is on the order of 1-100 years, the surface temperature abruptly plunges to 



values below 5 x 10^ K” (Lattimer et al. [167]). “Soon after a supernova explosion, the young neutron star has large 
temperature gradients in the inner part of the crust. While the powerful neutrino emission quickly cools the core, the 
crust stays hot. The heat gradually flows inward on a conduction time scale and the whole process can be thought of 
as a cooling wave propagation from the center toward the surface” (Gnedin et al. [170]). 

The two considerations we have quoted above are developed in the case of spherical symmetry and we would like to 
keep the mind open, in this new astrophysical field, to additional factors, some more traditional than others, to be taken 
into account. Among the traditional ones we recall: 1) the presence of rotation and magnetic field which may affect 
the thermal conductivity and the structure of the surface, as well as the above mentioned magnetospheric emission; 
2) there could be accretion of matter from the expanding nebula; and, among the nontraditional ones, we recall 3) 
some exciting theoretical possibilities advanced by Dyson on volcanoes on neutron stars [171] as well as iron helide 
on neutron star [172], as well as the possibility of piconuclear reactions on neutron star surface discussed in Lai & 
Salpeter [173]. 

All the above are just scientific arguments to attract attention on the abrupt fall in luminosity reported in this meeting 
on URCA-1 by Elena Plan which is therefore, in this light, of the greatest scientific interest and further analysis should 
be followed to check if a similar behavior will be found in future XMM and Chandra observations also in URCA-2. 


Astrophysical implications 

In addition to these very rich problematics in the field of theoretical physics and theoretical astrophysics, there are 
also more classical astronomical and astrophysical issues, which will need to be answered if indeed the observations 
of a young neutron star will be confirmed. An important issue to be addressed will be how the young neutron star can 
be observed, escaping from being buried under the expelled matter of the collapsing star. A possible explanation can 
originate from the binary nature of the newly born neutron star: the binary system being formed by the newly formed 
black hole and the triggered gravitational collapse of a companion evolved star leading, possibly, to a “kick” on and 
ejection of the newly born neutron star. Another possibility, also related to the binary nature of the system, is that the 
supernova progenitor star has been depleted of its outer layer by dynamic tidal effects. 

In addition, there are other topics in which our scenario can open new research directions in fundamental physics 
and astrophysics: 

1) The problem of the instability leading to the complete gravitational collapse of a ~ lOM© star needs the introduction 
of a new critical mass for gravitational collapse, which is quite different from the one for white dwarfs and neutron 
stars which has been widely discussed in the current literature (see e.g. Giacconi & Ruffini [6]). 

2) The issue of the trigger of the instability of gravitational collapse induced by the GRB on the progenitor star 
of the supernova or, vice versa, by the supernova on the progenitor star of the GRB needs accurate timing and the 
considerations of new relativistic phenomena. 

3) The general relativistic instability induced on a nearby star by the formation of a black hole needs some very basic 
new developments in the field of general relativity. 

Only a very preliminary work exists on this subject, by Jim Wilson and his collaborators, see e.g. the paper by 
Mathews and Wilson [174]. The reason for the complexity in answering such a question is simply stated: unlike the 
majority of theoretical work on black holes, which deals mainly with one-body solutions, we have to address here a 
many-body problem in general relativity. We are starting in these days to reconsider, in this framework, some classic 
work by Fermi [175], Hanni and Ruffini [176], Majumdar [177], Papapetrou [178], Parker et al. [179], Bini et al. [180] 
which may lead to a new understanding of general relativistic effects relevant to these astrophysical “triptychs”. 


The short GRBs 

After concluding the problematic of the long GRBs and their vast astrophysical implications, we have turned to 
the physics of short GRBs. We first report some progress in the understanding the dynamical phase of collapse, the 
mass-energy formula and the extraction of blackholic energy which have been motivated by the analysis of the short 
GRBs. In this context progress has also been accomplished on establishing an absolute lower limit to the irreducible 
mass of the black hole as well as on some critical considerations about the relations of general relativity and the 
second law of thermodynamics. We recall how this last issue has been one of the most debated in theoretical physics 
in the past thirty years due to the work of Bekenstein and Hawking. Following these conceptual progresses we analyze 



the vacuum polarization process around an overcritical collapsing shell. We evidence the existence of a separatrix 
and a dyadosphere trapping surface in the dynamics of the electron-positron plasma generated during the process 
of gravitational collapse. We then analyze, using recent progress in the solution of the Vlasov-Boltzmann-Maxwell 
system, the oscillation regime in the created electron-positron plasma and their rapid convergence to a thermalized 
spectrum. We conclude by making precise predictions for the spectra, the energy fluxes and characteristic time-scales 
of the radiation for short-bursts. 


Short GRBs as cosmological candles 

The characteristic spectra, time variabilities and luminosities of the electromagnetic signals from collapsing overcrit¬ 
ical stellar cores have been derived from first principles, and they agree with preliminary observations of short-bursts 
[19]. Hopefully new space missions will be planned, with temporal resolution down to fractions of /rs and higher 
collecting area and spectral resolution than at present, in order to verify the detailed agreement between our model and 
the observations. It is now clear that if our theoretical predictions will be confirmed, we will have a very powerful tool 
for cosmological observations: the independent information about luminosity, time-scale and spectrum can uniquely 
determine the mass, the electromagnetic structure and the distance from the observer of the collapsing core, see e.g. 
Fig. 43 and Ref. [54]. In that case short-bursts, in addition to give a detailed information on all general relativistic and 
relativistic field theory phenomena occurring in the approach to the horizon, may also become the best example of 
standard candles in cosmology [181]. We are currently analyzing the introduction of baryonic matter in the optically 
thick phase of the expansion of the plasma, within this detailed time-varying description of the gravitational 
collapse, which may affect the structure of the Proper-GRB (P-GRB) [13] as well as the structure of the long-bursts 
[28]. 


On the dyadosphere of Kerr-Newman hlack holes 

An interesting proposal was advanced in 2002 [182] that the plasma may have a fundamental role as well in the 
physical process generating jets in the extragalactic radio sources. The concept of dyadosphere originally introduced 
in Reissner-Nordstrom black hole in order to create the plasma relevant for GRBs can also be generalized to the 
process of vacuum polarization originating in a Kerr-Newman black hole due to magneto-hydrodynamical process of 
energy extraction (see e.g. [183] and references therein). The concept therefore introduced here becomes relevant for 
both the extraction of rotational and electromagnetic energy from the most general black hole [11]. 


H: H: * 

After the completion of these lectures, we have become aware that Ghirlanda et al. [184] have given evidence for 
the existence of an exponential cut off at high energies in the spectra of short bursts. We are currently comparing and 
contrasting these observational results with the predicted cut off in Fig. 42 which results from the existence of the 
separatrix introduced in [49]. The observational confirmation of the results presented in Fig. 42 would lead for the first 
time to the identification of a process of gravitational collapse and its general relativistic self-closure as seen from an 
asymptotic observer. 
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THE NUMERICAL INTEGRATION OF THE HYDRODYNAMICS AND THE RATE 
EQUATIONS IN THE LIVERMORE CODE 


The hydrodynamics and the rate equations for the plasma of e^e -pairs 

The evolution of the e+e^-pair plasma generated in the dyadosphere has been treated in two papers [55, 59]. We 
recall here the basic governing equations in the most general case in which the plasma fluid is composed of -pairs, 
photons and baryonic matter. The plasma is described by the stress-energy tensor 

+ (128) 

where p and p are respectively the total proper energy density and pressure in the comoving frame of the plasma fluid 
and is its four-velocity, satisfying 

8 niU‘f+8rriU^f = -l, (129) 

where U’’ and U‘ are the radial and temporal contravariant components of the 4-velocity and 

(f's = 8tt{r)d^t + 8rr{r)d^r -\-reef'd + r^siv^ dcf'(p , (130) 

where 8 tt{r) = —a^{r) and 8 rr{r) = 

The conservation law for baryon number can be expressed in terms of the proper baryon number density ng 

{nBU^)-g = 8^H8^nBU'')^v 

= {nBU‘)j + \{r^nBU''),r = 0. (131) 

The radial component of the energy-momentum conservation law of the plasma fluid reduces to 

^ + -l;((p + p)U'U,) + -^§-pr\p + p)U'U,)-'-{p + p)^^(U'f + ^{U'f'\ = 0 . ( 132 ) 


The component of the energy-momentum conservation law of the plasma fluid equation along a flow line is 

Ug{T^'').y = -{pU^).^-p{U^\^, 

= -g-kghu^),v-P8^kg^U'').v 

= (pt/0,) + 4(r2pt/'-),. 

ri 

+ p + =0. (133) 


We define also the total proper internal energy density e and the baryonic mass density Pb in the comoving frame 
of the plasma fluid, 

e = p-pB, Pb = ngme^ . (134) 


The numerical integration 


A computer code [185, 186] has been used to evolve the spherically symmetric general relativistic hydrodynamic 
equations starting from the dyadosphere [55]. 

We define the generalized gamma factor y and the radial 3-velocity in the laboratory frame V 
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From Eqs.(130, 129), we then have 


gtt a 
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Following Eq.(134), we also define 


(137) 


E = ey, D = PbJ, and p = py 

so that the conservation law of baryon number (131) can then be written as 


Eq.(133) then takes the form, 
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Defining the radial momentum density in the laboratory frame 

Sr=a{p + p)U‘Ur=iD + rE)Ur, 
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we can express the radial component of the energy-momentum conservation law given in Eq.(132) by 
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In order to determine the number-density of pairs, we use the pair rate equation. We define the 
pair density in the laboratory frame A(,± = yn^± and A^g±(r) = yng±{T), where n^±{T) is the total proper number 
density of pairs in comoving frame at thermodynamic equilibrium with temperature T in the process e+ -f ^ 7 -I- 7 
{ng-{m^T) =ny{T)'^, n^i is the total proper number density of pairs in comoving frame at a generic time before 
reaching the equilibrium. We write the rate equation in the form 


dNe± 

dt 


r^ dr' a 


N,±V^)+-^{N^±{T), 


(142) 


These equations are integrated starting from the dyadosphere distributions given in Eig. 17 (Right) in [28] and assuming 
as usual ingoing boundary conditions on the horizon of the black hole. 


THE NUMERICAL INTEGRATION OF THE HYDRODYNAMICS AND THE RATE 

EQUATIONS IN THE ROME CODE 

Era I: expansion of PEM-pulse 


After the explosion from the dyadosphere a thermal plasma of pairs and photons optically thick with respect 
to scattering processes begins to expand at ultrarelativistic velocity. In this era the expansion takes place in a region of 
very low baryonic contamination. 

Recalling that the local number density of electron and positron pairs created as a function of radius is given by 


ne+e-(r) 


Q 



(143) 



the limit on such baryonic contamination, where is the mass-energy density of baryons, is given by 


Pb, < nipiig+g- (r) = 3.2 • 




Near the horizon r ~ r+, this gives 


Pb, < mptig+g-^r) = 1.86 • 10 ^"^ 



(g/cm^), 


and near the radius of the dyadosphere 


PBc ^fnpng+g-{r) =3.2-10^ 1 



{g/crn^). 
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Such conditions can be easily satisfied in the collapse to a black hole, but not necessarily in a collapse to a neutron 
star. 

Consequently we have solved the equations governing a plasma composed solely of e+e^-pairs and electromagnetic 
radiation, starting at time zero from the dyadosphere configurations corresponding to constant density in Fig. 44. 

The plasma of pairs and photons is described by the covariant energy-momentum tensor: 




(147) 


where p and p are respectively total proper energy density and pressure in the comoving system; are contravariant 
components of 4-velocity and AT’^'' takes into account of dissipative effects due to heat conduction and viscosity, but 
in this treatment it has been neglected. In general we have gpvU^U'' = — 1. For a spherically symmetric motion this 
reduces to gtt{U‘)^ + grr{U'')'^ = — 1, where U‘ and t/'' are respectively temporal and radial controvariant components 
of 4-velocity U ^. 

It is assumed that the gravitational interaction with central black hole is negligible with respect to the total energy 
of PEM-pulse such that a fluid expansion with special relativistic equations can be considered. 

Moreover it is assumed that photons remain trapped inside fireball until complete transparency, i.e. the emission 
of electromagnetic radiation is negligible during the first phases of expansion, being therefore adiabatic [55]. This 
assumption is valid until the photon mean free path is negligible with respect to the thickness of pulse. 

The thermodynamic quantities used to describe the process are the total proper internal energy density of pulse e, 
given by e = 6 ^+ -I- e^- -I- Cy, where 6^+ (e^-) is total proper internal energy density of electrons (positrons) and Ey of 
photons. The proper number density of pairs if the system is in thermodynamic equilibrium initially at temperature 
T of order T ~ MeV, enough for pair creation, equals the proper number density of photons ny. This is not valid 
at lower temperature [156]. The pressure is p = p^+ -I- p^- + py, where p^± are electrons and positrons pressures and 
Py is photons pressure. The system is highly relativistic, so the equation of state ;? = e/3 can be considered valid. This 
equation of state is represented with thermal index F: 

r = i-f-. (148) 

e 


Fermi integrals 

Thermodynamical quantities before introduced are expressed in terms of integrals over Bose distribution for 
photons and Fermi distribution for pairs with zero chemical potentials Py and Pg±. We begin from the reaction 
^ Y + 7- From statistical mechanics it is known that given a thermodynamic system at temperature T kept 
inside a volume V and made of a number of particle variable N, the thermodynamic equilibrium is expressed by the 
condition that the potential free energy of Helmholtz F{T, V,N) is stationary with respect to N variations: 
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Figure 44. Three different dyadospheres corresponding to the same value and to different values of the two parameters jl and 
^ are given. The three different configurations are markedly different in their spatial extent as well as in their energy-density 
distribution (see text). 





















by definition chemical potential /x is given by 


(150) 


so that for a system made by a photon gas at equilibrium with matter with respect to creation and adsorption 
processes, we have jj.y = 0 [187]. Therefore the chemical potential of electrons and positrons associated to reaction 
^ 7 + yis equal and opposite: = —llg+ = fi; moreover also /r must be zero since the total electric charge 

of fireball is zero: if Q is total electric charge of fireball, we have 


Q = e[n^- {m,T,^i)-n^+ (m,r,-^)] =0 


where {m,T,jj.) is given by 
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In the following the expressions of thermodynamical quantities as Fermi integrals are listed. The proper number 
density of electrons [188] is given by 
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where z = pc/kT, m is the electron mass, ^[MeV] is the temperature of fireball in comoving frame, a is a constant 
given by a = 8 %^/\5h ^= 1.37 • Kp’^erg/cnr'MeV'^, k is the Boltzmann constant and A = (7/4)(7r^/15) is a 
numerical constant introduced for convenience. 

Since the thermodynamic equilibrium is assumed and in all cases considered the initial temperature is larger than 
pairs creation threshold (T = 1 MeV), the proper number density of electrons is roughly equal to that one of photons: 

~ (r) ~ Hy (T) ; (154) 

in these conditions the number of particles is conserved: 

K±t/'');^ = 0. (155) 

Later on, for T ^ IMeV (see Fig. 45), pairs go on in annihilation but can not be created anymore, therefore 

Hy (r) > n^± > Hgi (r) (156) 

as shown in Fig. 46. 

The total proper internal energy density for photons is given by 
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where p = hv/c. The total proper internal energy density for electrons is given by: 
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Figure 45. Temperature in comoving system as a function of emission time for different values of black hole mass fi. 


where z = pc/kT and the integral is computed numerically. Therefore the total proper internal energy density of the 
PEM-pulse, summing up all the contributions of photons and pairs, is given by 


£tot — tlT 


^ 7 2 + {mc^/kTy 

4 A Jo gs/z^+imc^/kTf _|_ 2 


(159) 


where the factor 2 in front of the integral takes into account of electrons and positrons. 
About the pressure of the photons it holds 



(160) 
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Figure 46. Ratio between number density of pairs Mg± and number density of photons ny (T) as a function of emission time 

for different values of black hole moss p. 


and about the pressure of electrons 
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Therefore the total pressure of PEM-pulse is given by 
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Numerical code 

In the following we recall a zeroth order approximation of the fully relativistic equations of the previous section 
[55]: 

(i) Since we are mainly interested in the expansion of the plasma away from the black hole, we neglect the 
gravitational interaction. 

(ii) We describe the expanding plasma by a special relativistic set of equations. 


In the PEM-pulse phase the expansion in vacuum is described by a set of equation expressing: 

• entropy conservation, because of the assumption that emission of electromagnetic radiation is negligible up to 
transparency; 

• energy conservation, because the increase of kinetic energy is compensated by a decrease of total internal energy. 


Eor the expansion of a single shell, the adiabaticity is given by 

d{Ve)+ pdV = dE + pdV = 0 


(163) 


where V is the volume of the shell in the comoving frame and E = Ve is the total proper internal energy of plasma. 
By using the equation of state 148 we find 

dlne + TdlnV = 0 (164) 

and, by integrating, we find 

e Vk 

recalling that the volume of the fireball in the comoving frame is given by V = where 'E is the volume in the 
laboratory frame, we find 


(165) 
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The total energy conservation of the shell implies [55]: 

(re)rf = (rso)ro^; 

and this gives the evolution for f: 


r=7o 


eo% 

ef 


Substituting this expression for 7 in (166) the final equation for proper internal energy density is found 


e = £0 


r 


(166) 


(167) 


(168) 


(169) 


The evolution of a plasma of pairs and photons should be treated by relativistic hydrodynamics equations 
describing the variation of the number of particles in the process. The 4-vector number density of pairs is defined 
(ng±t/^), which in the comoving frame reduces to the 4-vector (n^r, 0,0,0). The law of number conservation for pairs 
is 




1 






(170) 



where g =|| gf^'' ||= —r'^sin^O is the determinant of Reissner-Nordstr 0 m metric. In the system processes of creation 
and annihilation of particles occur due to collisions between particles. If the number of particles is conserved, it holds 
{ng±U^).^ = 0 ; if instead it is not conserved, in the assumptions that only binary collisions between particles occur 
and in the hypothesis of molecular caos, the Eq.(170) becomes 


{n,±U^).^ = av[n,-{T)n,+ {T) -n,-n^+] 


where a is the cross section for the process of creation and annihilation of pairs, given by 
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(172) 


with a = ^ and E total energy of positrons in the laboratory frame, and Ve 
the sound velocity in the fireball: 



the classical radius of electron, v is 

mc^ 


(173) 


and ^ is the mean value of crv; for a we use as a first approximation the Thomson cross section, Ot = 0.665 • 
lO^^^cm^; Hgi (T) is the total proper number density of electrons and positrons in comoving frame at thermodynamic 
equilibrium in the process —> 7 + 7 {n^- {m,T) = ny{T)), is the total proper number density of electrons 

and positrons in comoving frame at a generic time before reaching the equilibrium. 

Using the approximation of special relativity, the 4-velocity is written = (7,7^); substituting to n^± {T) the {T ) 

and to ng± the Eq.(171) in hybrid form becomes 
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dt 


lA 
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(174) 


valid for electrons and positrons. 

Now we have a complete set of equations for numerical integration: (169), (168) e la (174). 

If we now turn from a single shell to a finite distribution of shells, we can introduce the average values of the proper 
internal energy and pair number densities (e,ng±) for the PEM-pulse, where the average 7 -factor is defined by 


7 = 




(175) 


and 'E is the total volume of the shell in the laboratory frame [55]. 

In principle we could have an infinite number of possible schemes to define geometry of the expanding shell. Three 
different possible schemes have been proposed [55]: 


Sphere. An expansion with radial component of 4-velocity proportional to the distance to the black hole Ur{r) = 
r 

U , where U is the radial component of 4-velocity on the external surface of PEM-pulse (having radius 


^(f) 

the factor yfrom (175) is 
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this distribution corresponds to a uniform and time decreasing density, like in Eriedmann model for the universe; 

• Slab 1. An expansion with thickness of fireball constant ^ in laboratory frame in which the black hole 

is at rest, with Ur{r) =Ur = cost and j = a/ 1 -f Ur^; this distribution does not require an average; 

• Slab 2. An expansion with thickness of fireball constant in comoving frame of PEM-pulse. 

The result has been compared with the one of hydrodynamic equation in general relativity [55] (see Eig. 47). Excellent 
agreement has been found with the scheme in which the thickness of fireball is constant in laboratory frame: what 
happens is that the thickness in comoving frame increases, but due to the Lorentz contraction, it is kept constant 
in laboratory frame and equal to .^ = (r^s — ?"+)■ In this case Ur = \/ 7 ^-M, where y is computed by conservation 
equations. 

A similar situation occurs for the temperature of PEM-pulse. In the comoving frame the temperature decreases as 




Radius (cm) 


Figure 47. Lorentz y factor as a function of radial coordinate. Three schemes of expansion of PEM-pulse (see text) are compared 
with solution of hydrodynamics relativistic equations numerically integrated for a black hole with fi = 10^ and ^ =0.1. The result 
is in accordance with the scheme of a fireball with constant thickness in laboratory frame. 


T' in accordance with results in literature [29]. Since 7 monotonically increases as y ^ R [13], in laboratory 

frame T — jT' ~ constant [59]; photons are blue-shifted in laboratory frame in such away that, at least in the first 
phase, the temperature measured by an observer at infinity is constant. The numerical value of the temperature of 
equilibrium at each instant is found by imposing the equivalence, within a certain precision, of (159) numerically 
computed and (169). 

Even if the PEM-pulse is optically thick in the expansion before transparency, photons located at a distance from 
the external surface less their mean free path can escape and reach the observer at infinity. The mean free path in the 
comoving frame is given by 

Ly =—-lO^'^cm (177) 

while in laboratory frame is given by A = Lylf ^ lO^^cm. However the luminosity emitted at this stage is negligible, 
since the ratio between X and the thickness of the fireball & in the laboratory frame (with & = [r^s — r+) ~ 10 ®cm) is 




of the order of ~ 10 


Era II: interaction of the PEM pulse with remnant 


The PEM pulse expands initially in a region of very low baryonic contamination created by the process of 
gravitational collapse. As it moves outside the baryonic remnant of the progenitor star is swept up. The existence 
of such a remnant is necessary in order to guarantee the overall charge neutrality of the system: the collapsing core 
has the opposite charge of the remnant and the system as a whole is clearly neutral. The number of extra charges in 
the baryonic remnant negligibly affects the overall charge neutrality of the PEM pulse. 

The baryonic matter remnant is assumed to be distributed well outside the dyadosphere in a shell of thickness A 
between an inner radius and an outer radius rout = rjn + A at a distance from the black hole not so big that the PEM 
pulse expanding in vacuum has not yet reached transparency and not so small that the system will reach enoughly high 
value of Lorentz 7 in order to not be stopped in the collision (see Eig. 9). Eor the sake of an example we choose 

Tin = lOOrds, A = lOrds- (178) 


The total baryonic mass Mb — Ngrrip is assumed to be a fraction of the dyadosphere initial total energy (Edya)- The 
total baryon-number Ab is then expressed as a function of the dimensionless parameter B given by 


NBrripC^ 

Ttdya 


(179) 


where B is a parameter in the range 10^^ — 10^^ and mp is the proton mass. We shall see below the role of B in the 
determination of the features of the GRBs. We saw in section the sense in which B and can be considered to be 
the only two free parameters of the black hole theory for the entire GRB family, the so called “long bursts”. Eor the 
so called “short bursts” the black hole theory depends on the two other parameters /r, ^, since in that case B — 0 since 
most of the energy, unless the whole energy, in the pulse is emitted at transparency. The baryon number density tip is 
assumed to be a constant 

—o —o —o 2 /I on\ 

= PB = mpnBC . (180) 

Vb 

As the PEM pulse reaches the region < r < rout, it interacts with the baryonic matter which is assumed to be at 
rest. In our model we make the following assumptions to describe this interaction: 

• the PEM pulse does not change its geometry during the interaction; 

• the collision between the PEM pulse and the baryonic matter is assumed to be inelastic, 

• the baryonic matter reaches thermal equilibrium with the photons and pairs of the PEM pulse. 

These assumptions are valid if: (i) the total energy of the PEM pulse is much larger than the total mass-energy of 
baryonic matter Mb, 10^^ < B < 10^^, (ii) the ratio of the comoving number density of pairs and baryons at the 
moment of collision n^+g- /rig is very high (e.g., 10® < n^+g- /n°g < 10*^) and (iii) the PEM pulse has a large value of 
the gamma factor (100 < 7 ). 

In the collision between the PEM pulse and the baryonic matter at rout > ?" > ^'in , we impose total conservation of 
energy and momentum. We consider the collision process between two radii r 2 ,ri satisfying rout > r 2 > ri > rtn and 
r 2 — ri ^ A. The amount of baryonic mass acquired by the PEM pulse is 


MbTtt 3 

AM =-(r, — r 

Eb 3 ^ 2 


( 181 ) 


where Mb/Vb is the mean-density of baryonic matter at rest. 

As for energy density of dyadosphere, here also we choose a simplification for the energy density: in fact during the 
passage of the shell a deposition of material on the external surface of the fireball creates; however we neglected this 
effect and assumed that this material after collision diffuses instantaneously in the pulse with a constant density: 


( 182 ) 



where N'^ is the number of particle of the remnant shell swept up by the pulse and V is the comoving volume of the 
fireball. 

The conservation of total energy leads to the estimate of the corresponding quantities before (with “o”) and after 
such a collision 

(TEo + p^)^%+AM= (Te + Pb+ ^+rAs)fr, (183) 

where Ae is the corresponding increase of internal energy due to the collision. Similarly the momentum-conservation 
gives 

(TEo + p°B)foU°% = {re + pB + ^ + rA£) 7 t/,r, ( 184 ) 

where the radial component of the four-velocity of the PEM pulse is U° = 1 / 7 ^ — 1 and E is the thermal index. We 
then find 
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(185) 

(186) 


These equations determine the gamma factor 7 and the internal energy density £ = £0 + Ae in the capture process of 
baryonic matter by the PEM pulse. 

The effect of the collision of the PEM pulse with the remnant leads to the following consequences: 

• a reheating of the plasma in the comoving frame but not in the laboratory frame; an increase of the number of 

pairs and of free electrons originated from the ionization of those atoms remained in the baryonic remnant; 
correspondingly this gives an overall increase of the opacity of the pulse; 

• the more the amount of baryonic matter swept up, the more internal energy of the PEMB pulse is converted in 
kinetic energy of baryons. 

By describing the interaction of PEM pulse with remnant as completely inelastic collision of two particles, one can 
compute by the energy-momentum conservation equation the decrease of Lorentz 7 and the increase of internal energy 
as function of B parameter and also the ultrarelativistic approximation ( 7 , ^ 

1 . an abrupt decrease of the gamma factor given by 

_ 1+B B+1 

■ ’'Vr.-(2fi+fi2) + i TWW • 

where 70 is the gamma factor of the PEM pulse before the collision, 

2 . an increase of the internal energy in the comoving frame £^ 0 // developed in the collision given by 


Ecoll y/7o^(2B+B2) + l 

^dya To 


— >y^^oo-B+VB^ + 2B, 


This approximation applies when the final gamma factor at the end of the PEM pulse era is larger than Jcoih right panel 
in Eig. 9. 

In this phase of expansion, another thermodynamic quantity has not been considered: the chemical potential p of 
the electrons from ionization of baryonic remnant. We remind that the total proper number density of electrons of 
ionization is given by 
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four equations are imposed to find a formula useful for numerical computation: the first one is the thermodynamical 
equilibrium of fireball, or 

«e±(7’o) = 


Hg- = 'Z.riB 


the second one is 


(189) 



where 1/2 < Z < 1, with Z = 1 for hydrogen atoms and Z = 1/2 for baryonic matter in general; the third one derives 
from the definition of B, and states a relation between the two densities hs and from definition of B, we have 


Nb 
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(190) 


where To is the initial temperature of fireball and is a parameter (b < 0) defined by (190); so if Vo is the initial volume 
of dyadosphere and w the initial volume of the baryonic shell 


4 = 10*n,±(ro) —; 

W 

finally the fourth one is the conservation law of baryonic matter 

(n^t/^);^=0. 


(191) 


(192) 


Therefore the chemical potential /r is numerically determined at a certain time of expansion if the initial temperature 
To of fireball and the initial volume of baryonic shell w are known and, at that time, the volume V, the temperature T 

n‘’_{T) 

and the Lorentz factor 7 of the fireball, the volume of the baryonic shell swept up vb and the ratio —: 




^y/z^ + {mc^/kT)^+-^ _|_ j 


dz 


(193) 


where the factor in brackets must be considered only for r > rout, while the proportionality factor is the function 
zeta of Riemann 1 ^(x) for computation of ny, with i^(3) = 1.202. 

Therefore the equations for this phase are (185), (186), (182), (174) and (193). 


Era III: expansion of PEMB pulse 


After the engulfment of the baryonic matter of the remnant the plasma formed of -pairs, electromagnetic 
radiation and baryonic matter expands again as a sharp pulse, namely the PEMB pulse. The calculation is continued 
as the plasma fluid expands, cools and the pairs recombine until it becomes optically thin; 


L 


dr{n^± +ZnB)GT — 0(1), 


(194) 


where Gt = 0.665 • lO^^^cm^ is the Thomson cross-section and the integration is over the radial interval of the PEMB 
pulse in the comoving frame. In order to study the PEMB pulse expansion the validity of the slab approximation 
adopted for the PEM pulse phase has to be verified; otherwise the full hydrodynamics relativistic equations should be 
integrated. The PEMB pulse evolution firstly has been simulated by integrating the general relativistic hydrodynamical 
equations with the Livermore codes, for a total energy in the dyadosphere of 3.1 x 10^“^ erg and a baryonic shell of 
thickness A = lOrds at rest at a radius of lOOrds and B ~ 1.3-10^^. 

In analogy with the special relativistic treatment for the PEM pulse, presented in section (see also [55]), for the 
adiabatic expansion of the PEMB pulse in the constant-slab approximation described by the Rome codes the following 
hydrodynamical equations with ps ^0 has been found 
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Figure 48. Lorentz y factor as a function of radial coordinate from the PEMB-pulse simulation is compared with the y factor as 
solution of hydrodynamics relativistic equations numerically integrated (open squares) for E^ya = 3.1 x lO^^erg and 6 = 1.3 x 10^^, 
rjn = lOOrjs and A = lOrjs. The result is in accordance with the scheme of a fireball with constant thickness in laboratory frame 
which is valid up to B = 10^^. 


In these equations (r > rout) the comoving baryonic mass and number densities are ps — Mb/V and ng = Nb/V, where 
V is the comoving volume of the PEMB pulse. 

The result is shown in Fig. 48 [59] where the bulk gamma factor as computed from the Rome and Livermore codes 
are compared and very good agreement has been found. This validates the constant-thickness approximation in the 
case of the PEMB pulse as well. On this basis we easily estimate a variety of physical quantities for an entire range of 
values of B. 

For the same black hole different cases have been considered [59]. The results of the integration show that for the 
first parameter range the PEMB pulse propagates as a sharp pulse of constant thickness in the laboratory frame, but 
already for B ~ 1.3-10^^ the expansion of the PEMB pulse becomes much more complex, turbulence phenomena can 
not be neglected any more and the constant-thickness approximation ceases to be valid. 

It is also interesting to evaluate the final value of the gamma factor of the PEMB pulse when the transparency 
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Figure 49. Left) The gamma factors are given as functions of the radius in units of the dyadosphere radius for selected values 
of B for the typical case E^ya = 3.1 x 10^“* erg. The asymptotic values 7asym = £dya/(^BC^) = 10“*, 10^, 10^ are also plotted. The 
collision of the PEM pulse with the baryonic remnant occurs at r/r^^ =100 where the jump occMrj. Right) The y factor (the solid 
line) at the transparency point is plotted as a function of the B parameter. The asymptotic value (the dashed line) E^y^/(Mgc^) is 
also plotted. 


condition given by Eq.(194) is reached as a function of B, see Fig. 49. For a given black hole, there is a maximum 
value of the gamma factor at transparency. By further increasing the value of B the entire E^ya is transferred into the 
kinetic energy of the baryons (see also [59]). 

In Fig. 49Feft we plot the gamma factor of the PFMB pulse as a function of radial distance for different amounts of 
baryonic matter. The diagram extends to values of the radial coordinate at which the transparency condition given by 
Fq.(194) is reached. The “asymptotic” gamma factor 
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is also shown for each curve. The closer the gamma value approaches the “asymptotic” value (199) at transparency, 
the smaller the intensity of the radiation emitted in the burst and the larger the amount of kinetic energy left in the 
baryonic matter (see Fig. 49Right). 


ON URCA PROCESS 


From G. Gamow [189]: 

“The summer of 1939 I spent with my family vacation on the Copacabana beach in Rio de Janeiro. One evening, 
visiting the famous Casino da Urea to watch the gamblers, I was introduced to a young theoretical physicist born on 
an Amazon River plantation, named Mario Schoenberg. We became friends, and I arranged for him a Guggenheim 
fellowship to spend a year in Washington to work with me in nuclear astrophysics. His visit was very successful, and 
we hit upon a process which could be responsible for the vast stellar explosions known as supernovae. The trick is done 
by alternative absorption and reemission of one of the thermal electrons in the very hot (billions of degrees!) stellar 
interior by various atomic nucleai. Both processes are accompanied by the emission of neutrinos and antineutrinos 
which, possessing tremendous penetrating power, pass through the body of a star like a swarm of musquitoes through 
chicken wire and carry with them large amount of energy. Thus, the stellar interior cools rapidly, the pressure drops, 
and the stellar body collapse with a great explosion of light and heat. 

All this is too complicate to explain in non technical words, and I am mentioning it only as background for how we 
came to give that process is name. We called it the Urea process, partially to commemorate the casino in which we first 










met, and partially because the Urea process results in a rapid disappearance of thermal energy from the interior of the 
star, similar to the rapid disappearance of money from the pockets of the gamblers of the Casino d Urea. Sending our 
article “On the Urea process” for publication on the physical review I was worried that the Editor would ask why we 
called the process “Urea”. After much thought I decided to say that this is short for “UnRecordable Cooling Agent”, 
but they never asked. Today, there are other known cooling processes involving neutrinos which work even faster than 
the Urea process. For example, a neutrino pair can be formes instead of two gamma quanta in the annihilation of a 
positive and negative electrons”. 


CASINO DE URCA TODAY 
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